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ABSTRACT 



A proposed neutron**neutron scattering experiment in which the neutron 
source would be colliding beams from an underground nuclear explosion^was 
investigated. In particular, the accuracy to be expected was determined 
so as to allow evaluation of the worth of the experiment. Phase shift 
analysis was performed for Gaussian, Exponential, Yukawa, and Square 
Weil spherically symmetric potentials. The energy range considered was 
20 kev to 2.6 MeV. Scattering length was found to be determined to 
within 6% and effective range within k0% for a 10% scattering cross section 
accuracy. Potential depth would be determined to at least .5% and range to 
.3%. 
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1. Introduction 



A low energy neutron-neutron scattering experiment utilizing a 
single underground nuclear explosion has been proposed and studied. 
Preliminary studies indicate the scattering cross section might be meas- 
ured to an accuracy of 10% from 20 kev to 2HeV. The efforts of this 
paper are directed to a numerical investigation of the expected scatter- 
ing with special attention given to: (a) the possible existence of the 

dineutron; (b) the effect of changes in effective range and scattering 
length upon the range and depth of the spherically symmetric Gaussian, 
Exponential, Yukawa, and Square Well Potentials; (c) the effect of changes 
in effective range and scattering length upon the scattering cross 
section from each of the potential shapes in (b) ; and (d) a comparison 
of calculated cross sections from the four potential shapes in (b) with 
the shape-independent cross section. 

Precise determination of the neutron-neutron scattering parameters 
would be important for a number of reasons as stated by Moravcsik. Q]] 
These reasons Include: (1) an investigation of charge independence; (2) 

an investigation of charge symmetry; and (3) the existence or non-existence 
of a slightly unbound dineutron resonance state. 

This investigation attempted to locate a dlneutron state in section 
2 to provide a theoretical justification for its experimental quest. In 
Sections 3 and 4, 2 to 16% deviations in the ratio of effective range to 
scattering length from the accepted n-p ratio were investigated for their 
effect on scattering cross sections and Gaussian, Exponential, Yukawa, 
and Square Well Potential parameters. These results should allow some 
determination of the accuracy in experimental results required to 
determine to a given accuracy the shape-independent scattering parameters. 
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In Section 5, the cross section from effective-range theory is compared 
with cross sections from potential scattering to test the validity of 
the shape independent approximation. This report concludes with a 
summary in section 6 and the complete computer program used in the 
calculations along with the calculation results as a part of the appendix. 

The analysis was performed using phase shift scattering methods. 

The well established scattering length and effective range for neutron- 
proton scattering, a*^ and r^ were taken as a base for calculations. 
The ratio of effective range to scattering length was formed to yield 
a spherically symmetric potential range and depth for square well, 
exponential, gaussian, and yukawa potentials. The absolute value of 
the scattering length was then varied in 2% steps up to l6% while 
holding the effective range constant at its accepted value. 17 sets 
of potential parameters were thus obtained for each of the four poten- 
tial shapes. Then, using the maximum and minimum acceptable value of 
effective range as a base, the effective range was varied in 2% Incre- 
ments up to 6% while holding the scattering length constant at its 
accepted value. Eight additional sets of potential parameters were 
obtained for each of the potential shapes. 

The spherically symmetric radial wave equation was then numerically 
integrated using a fourth order runge kutta method for all calculated 
sets of parameters and potential shapes. Integration was performed for 
30 values of center of mass energy ranging from 20 kev to 2.6MeV and 
was carried out to a radius at which the effect of the nuclear potential 
was no longer felt. 

The logrithmic derivative of the numerically integrated solution 
was then compared to the logrithmic derivative of the exact field-free 
asymptotic solution outside the range of the nuclear potential and a 
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phase shift obtained. 

The total cross section as a function of center of mass energy was 
then calculated from the phase shift for each case. A ratio of each 
cross section to that obtained from the accepted neutron-proton scatt- 
ering potential parameters was formed. This gave the per cent change 
in total cross section for each incremental change in scattering length 
and effective range. 

Per cent change in total cross section from that for n-p scattering 
parameter calculations was plotted against energy for significa<nt incre- 
ments In two different energy ranges, 20 kev to 100 kev and 100 kev to 
2.6MeV. Potential depth and range were plotted against the ratio of 
effective range to scattering length for each potential shape. A com- 
parison of shape- independent effective range scattering calculations 
with potential scattering calculations was made by plotting the total 
cross section for each case on the same graph. 



3 




0 



2. Possible dineutron state 



An attempt was made to justify the possible existence of a dineutron 
resonance state near zero energy using a compound nucleus model, but 
failure was met In attempting to formulate a compound nucleus model for 
neutron-neutron scattering. 

Numerical calculations for neutron-neutron potential scattering 
were made for a large number of potential depths on the Control Data 
Corporation I604 digital computer in an attempt to find possible 
resonances, particularly a sharp resonance. Gaussian, exponential, 
and yukawa potentials were used with potential depths ranging from 
0 to the value obtained from the accepted n-p scattering parameters 
in 2.5 MeV steps over two center of mass energy ranges, 20 to 100 kev 
and .1 to 2.6 MeV. Results were scanned for resonances. No reson- 
ances were expected and none were found. 
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3. Effective range and potential parameters 



Detailed comparisons with experiments have given for the singlet 
effective range and scattering length for neutron-proton scattering: 



For an examination of the effect of deviations from this scattering 
length on the scattering cross section, the ratio ra,^ formed 

and a^ was varied in magnitude from 0-16% In 2% steps while holding r^^^ 
at its accepted value, 2.46 Fm. To examine the effect of deviations in 
effective range on the scattering cross section, the ratio r^^g / 
was again formed using the accepted value for a^, -23.S77 Fm, and the 
maximum and minimum values for r^^ ,2.58 Fm and 2.34 Fm respectively. 

**0,8 varied up to 6% in 2% steps in the direction of increasing 
deviation from the accepted value, 2.46 Fm. 

Four spherically symmetric potential shapes were then examined. 
These potentials are: 



6± 0.12 FM = - 23 . 6 77 ± 0.02<l FM L4Q 



Square Well: V = V© Fn 



V - 0 



r>Fn 



Gaussian Wei I : 



V = -V^ 






Yukawa Wei I : 



- Vc 




Exponential Wei I : 



V = - Vo C' F/rnl 
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Blatt and Jackson have arrived at a correspondence between the 
ratio of effective range to scattering length and the potential para- 
meters, Vq and r^, for each of the well shapes investigated. QQ 
The well parameters are given in terms of the intrinsic range, b, and 
the well depth parameter, S. The well depth, for the potential 

in question Is multiplied by a constant such that the first resonance 
occurs at zero energy. The well depth parameter Is then the inverse 
of this constant given by: 

Vcr)= 

The intrinsic range, b, is determined by the shape and range of the 
adjusted potential through the radial wave function for that potential. 
I t is defined by: 

- ulcr)3 dr 

kES o 

Blatt and Jackson provide interpolation formulas to relate and 
r„ to r^ and a. The form of their formulas have been changed slightly 
to fit the definition of the potentials in this report. In their 
modified form, these Interpolation formulas are: 

Square Wei 1 : 

s a C.I- 1 . 574<r,/a) v.35aCr*/o.>^3/D“ i- 1.3 0 (ro/cO**] 

Vo = l.0^^ 76 )O0 S / 

r„ ^ b 

Gaussian Wei 1 : 

6 5 LI “2, 030 Cn,/(X) + .567 (re/ 0 L>^ 3 / LI - £. 8»0 Cn./0-)4l,©05 C.ro/ci >^2 
b =r6/Cl-.64Us-0 C6-l?4U7t CS-l>’‘-i.Ol5-(S-|f J 

Vo*- aos X S/b^ 
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r„ = t b* /£. 0 6O4l''^ 

Exponential Well: 

S = CI- + 1.3 6oa (r«>/cx.f □ /Cl- 3.H34tro/0.') + 'i.u T7(n>/a)J 

1}= r*/LI-.3o4C5.-0 +.145<.s-»^^-.&43C^“0®+.n4 

Vo- 7.SIS-4I 

r„ = t/3.s-4ia 

Yukawa Well: 

SiC»- 3.ffoq3<n>/OL>+3.o 35‘4<'ro/<jj?' j/ci- » 34 3 cr<>/0w> + ^.loiaCro/a’Jj 
w» rj/c.»- 1. (a -n + 1 . 0 33(6-0^- \ . \ 8.nC3-i')® 1 1 . oorc%- 1'?’-. 3SL4(3-\f 3 
Ve = 3.l^9^ U66 a /b** 

fn «• b/2.,||Cife 

The results of these calculations are presented both In tabular and 
in graphical form In tables I and II and In graphs I through l6. The 
tables are arranged to exhibit and r^| as a function of the percentage 
deviation of r^ and a from r*^ ^ and a*^, the best values from n-p 
scattering. The graphs depict Vq and r,, as a function of r^/^ using 
the left and bottom scales respectively. Using the top and right 
scales, the percentage deviation of Vq and r„ from the best n-p scatter- 
ing parameters is shown as a function of the percentage deviation of 
ro/a from 

Variations In potential parameters were found to be of the order 
of a .5% increase In IVqI and a .5% decrease In r^ for a 10% increase 
In scattering length. A 10% increase In effective range, however, was 
observed to cause a decrease In IVqI of approximately 8% while Increasing 
r^ about 3%. This data is shown in table III. 
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EFFECT OF CHANGES IN EFFECTIVE RANGE UPON POTENTIAL PARAMETERS 
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TABLE III 



PERCENTAGE CHANGE IN WELL PARAMETERS CAUSED BY A TEN PERCENT INCREASE 
IN SCATTERING LENGTH AND EFFECTIVE RANGE 
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4. Phase shift analysis 



A partial wave analysis of neutron-neutron scattering was performed 
for ail potential parameters calculated in section 3« 

Assuming a spherically symmetric potential, the radial wave equa- 
tion 

R- O (I) 

becomes 

4- U ^ EUL (2) 

where U(r^^rRCr> 

Rearranging, equation (2) becomes 

U (3) 

where H E = C E ( m6V> / 4 ^ I 

and V (r> = H V(r) X 10^^ Q6 ] 

Ms neutron-proton reduced mass 
E= center of mass energy z 1/2 laboratory energy 
Separating equation (3) Into two first order equations for simul- 
taneous numerical integration gives 

^ -- - V CO - a <5) 



27 



Equations (4) and (5) were numerically integrated using a fourth 



order Runge-Kutta method with Boundry conditions: 

For r = O 

UCO) -O 

Z (O) = ± 

integration was carried outside the range of the nuclear potentials 
with 500 to 750 steps to about 15 fm for two center of mass energy 
ranges; 0-100 kev in 20 kev intervals, and .1-2.6 MeV in .1 MeV intervals. 
The values of Z and U at r i 15 fm are referred to as Zp and Up. 

The exact asymptotic solution of the field-free wave equation, 
equation (1) with V ■ 0, is given by Zeleny as: 

RjCr>=A 5 ^cos jjUr), - sin r\jj CSkrtJ (6) 

from which 






With the substitutions 



a, - r 




equations (6) and (7) become: 



U^(0= A^r{cos s, cgtrt - ^\n 



( 8 ) 



= A| (rC ^ CO s Sjj (J - (HO 

d r 



( 9 ) 



sin (I n^_,(hO-CH0n^^j(.»?r))} 
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where Is the scattering phase shift and and n^ are spherical 
Bessel and Neumman functions respectively. 

The scattering phase shift, ^ , Is obtained by matching the 
logrithmic deviations of the exact field-free and numerical solutions 
at a range where the nuclear force Is negligible. This procedure 
yields ; 

For : 0 

•= n^CKr) cos Sft 

r 

r»»s PM 

c^r) sin CoS <5* - riot<fep^ sin 

For X = ] 
r»15FW 

+ cosS,i,(0H*>-smS,n,(«tr5jArtosft, j,<l?r)-S\n ^C^^rSS\l 

»r'i5 PP| 

Rearranging and taking the Inverse tangent we get the following express- 
ions for the phase shift: 

- -toLn"' 



jKn,(kf^]r - HoC^r) J 



( 10 ) 



§, = t an' < 



n, Cfer)}r +niCfer^ 



( 11 ) 



The total cross section as a function of energy Is obtained from 
the phase shift using the relation 



=5-1^ C?.Q+0 ^trv^ So 



( 12 ) 
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For Identical particle scattering, equation ( 1 21 ) becomes, for the 



singlet state: 

1 ^ ^ I 

even ^ 

and for the triplet state: 




The computation of cross sections was made for both the singlet 
^ s 0 state and the triplet = 1 state. The : 1 cross sections 
were observed to be at least 3 orders of magnitude smaller In the 
energy range under consideration and were disregarded In subsequent 
computations. The computer program In appendix Ml however, retains 
the capability of ■ 1 calculations. 

The computer program In appendix Ml was used In cross section 
calculations for all the potential parameters In section 3. First 
the program computed the potential parameters at the base scattering 
length and effective range and at percentage deviations from this base 
read Into It In Its Input section. The total scattering cross section 
was then computed as a function of center of mass energy. The calcul- 
ated cross section for each set of potential parameters was then divided 
by the cross section for the base parameters to normalize the cross 
section and present It as a deviation from the base. 

The resulting data was plotted to show the percentage difference 
in total cross section from the calculated total cross section for n-p 
scattering with a^ ■ - 23.677 fm and s 2.46 fm as a function of 
energy for percentage deviations of a and r^. The graphs are drawn 
over two distinct energy ranges, 20 - 100 kev and .1-2.6 MeV, to 
better portray results and separate graphs are provided for each 
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potential shape. 

In Moravcsik's paper, the function Q(T) is defined by: 

where g » MT for T ■ laboratory kinetic energy and m s neutron mass. 

2 . 

Horavcsik has determined that in a 10% experiment, a linear approxi- 
mation for Q(T) would be sufficient. 

The value of this straight line at £ s 0 would give 1/^2 directly 
while the slope of the straight line would determine and the sign 

of a. 

Moravcsik concludes that a determination of total cross section 
to 10% accuracy at ten suitably distributed energy points could det- 
ermine the scattering length to 3% and the effective range to 50 - 70%. 

Tables IV, V of this report show that a 10% determination of 
total cross section in the vicinity of the n-p scattering parameters 
will allow a 6% determination of a$ at 20 kev and a 16% determination 
at 400 kev while at 2.6 MeV a 16% error in scattering length has only 
a 2% effect on total cross section. The same tables show that a 10% 
determination of total cross section will allow a 40% determination 
in effective range at 2.6MeV while a 100% error in effective range 
has approximately a 2% effect at 20 kev. 

An alternate method of determining the effective range parameters 
is therefore offered. The magnitude of the total cross section at 
20 kev could be compared to supplemented data of this report to arrive 
at a value for a^ within 6%. Graph 17 is shown as an illustration of 
possible supplemental data. The magnitude of the total cross section 
at 2MeV or greater could, in the same manner, be used to arrive at a 
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TABLE IV 



EFFECT OF A 10% ERROR IN TOTAL CROSS SECTION UPON THE ERROR IN a, 

9 



% Error In for 



ENERGY (kev) 


SQUARE WELL 


GAUSS IAN 


EXPONENTIAL 


YUKAWA 


20 


8-16 


6 


6 


6.5 


60 


8-16 


8.5 


8.5 


9.5 


100 


8-16 


11 


12 


16 


200 


8-16 


16.5 


16.5 


18 
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TABLE V 



EFFECT OF A 10% ERROR IN TOTAL CROSS SECTIONS UPON THE ERROR IN 

%error In r for 



ENERGY (MeV) 


GAUSS IAN 


EXPONENTIAL 


YUKAWA 


1.6 


60 


58 


62 


1.8 


56 


52 


59 


2.0 


53 


50 


57 


2.2 


50 


47 


55 


2.4 


45 


45 


52 


2.6 


40 


42 


50 
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value for r_ - within 50%. Because of the negligible effect of both 
large errors in r^^^ at the low energy cross section and large errors 
in a^ at the high energy cross section, they may be determined indepen* 
dently, each having little effect upon the determination of the other. 
Either gaussian, exponential, or yukawa data could be used because of 
only about a 1% difference in a^ value among them for any given cross 
section. Thus the scattering length and effective range could be det- 
ermined with a minimum of two points, one at either end of the energy 
range. For greater accuracy, the computer program in the appendix 
could easily be adapted to a searching interpolation. An advantage 
could be gained if the necessity for only two values of cross section 
would allow a more precise determination of those two values. It can 
be seen from accompanying graphs that the energy must be accurate to 
10% to achieve meaningful results in this way. If this energy accuracy 
is not feasible, then the method outlined by Horavcsik for scattering 
length and effective range determination must be used. In any event, 
Horavcsik's development of expected accuracies is confirmed by this 
phase shift analysis. 



35 



i* 



• b<fPET^ENCE. W Total O^qs^ 




GRAPH 18 • 





37 



f 





m 




GRAPH 20 

38 



t 



I N ~roTAU,_CR .QS.S--5iLCrJ-£>Ai.I^^ 






40 



1 . T Q_TA L ■ ■ f^Q S5: <^e C.T t OM 




t (KEV) 



GRAPH 23 

IjI 



J3JEE££aiic£,„wjjciTya Section C ^) 












43 







GRAPH 26 



44 




GRAPH 27 











I 







15 



I 

i 




GRAPH 28 



46 



PJF-E■BR£^j.C.P.^_l>J TOT/^ Cl?Q^S 




47 





# 



o 




GRAPH 30 



48 




J)i£EE-R£NJC6 tjJJToT Al— C l?QSS ^E CTloM ( 




GRAPH 31 



S. Shape independent approximation 



The total cross section for neutron-neutron scattering was cal- 
culated by the shape independent approximation of effective range 
scattering theory using r*^ ^ and This cross section was compared 

with that from the potential scattering calculations in section 4. 



where a(k) is the scattering length as a function of wave number. An 
approximation for a(k) is given for k ^ 0 as: 



which, when substituted for a(k) in the cross section formula yields: 



Formula (la) was used in a calculation of total cross section in 
the energy range .1 - 2.6MeV and compared with the cross sections 
obtained from potential scattering for the same parameters. The results 
are shown in graphs 32 - 33. 

The results of the effective range theory calculations compare 
favorably with potential scattering calculations. A 1% difference is 
the greatest deviation noted In the comparison. This confirms the 
commonly held belief that the shape independent approximation of 
effective range theory is quite accurate in this energy range. 



/ 



Segre gives for the total cross section 
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6. Sunmary 



In section 4 ft was determined that a 10% scattering experiment 
could define a scattering length to approximately 6% and an effective 
range to approximately 40% accuracy. This accuracy might be achieved 
with only two values of cross section; that at 20 kev and that at 2.6 
HeV, If the energy could be determined to 10% In both regions. Certainly 
It would be achieved by a least squares fit at 10 - 20 appropriately 
distributed energies as proposed by Horavcsik. 

In the region of the accepted n-p scattering length and effective 
range, It was found that a 2% change In scattering length resulted In 
a 3«1% change In cross section at 20 kev center of mass energy while 
only causing a .1% change In total cross section at 2.6MeV. A 20% 
change In effective range was seen to cause a 5% change In cross section 
at 2.6MeV while having almost no effect at 20 kev. Thus the conclusion 
was reached that effective range and scattering length could be determined 
separately and Independently by the cross sections at the high and low 
ends of the proposed energy range. 

In section 3 it was found that a 1% change In potential depth 
or range resulted in about a 20% change in scattering length and about 
a 1% change In effective range for Exponential and Gaussian wells. 

For the Yukawa well, a 1% change In potential depth and range 
resulted in about a 12% change in scattering length and a 1% change 
In effective range. 

It was conversely determined that a 6% determination of scattering 
length can determine the Gaussian and Exponential potential depths to 
less than - .2% and the potential range to less than - .3% for the 

n-p scattering length and effective range vicinity. The Yukawa well 
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depth would be determined to less than -.5% and the potential range 
to -.2%. 

It therefore seems quite likely that a 10% scattering experiment 
would determine the scattering length and effective range accurately 
enough to furnish an acceptable set of potential parameters. These 
potential parameters appear to be adequate to investigate charge - 
independent and charge-symmetry effects. 

Effective range theory was quite adequate to determine cross 
sections in the examined energy range. A 1% agreement was found with 
potential scattering. There should be little question of its propriety 
here. 
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APPENDIX I 



COMPUTER PROGRAM DESCRIPTION 

The program Is written In FORTRAN II language and, as listed In 
Appendix III, Is arranged to operate with the Control Data Corporation- 
1604 as a central processor and the IBM- 140I as peripheral equipment. 
The library tape functions used are the standard FORTRAN II functions 
with the exception of a subroutine called DRAW. This subroutine Is a 
part of the library tape at the Computer Center of the United States 
Naval Postgraduate School and is identified by Its large computer 
operators co-op ID; J7“NPS-DRAW. The subroutine was written by J.R. 

Ward of the U.S. Naval Postgraduate School. Where this subroutine Is 
not available, the program can be run by omitting subroutines SKETCH 
and NOW and their call cards. 

The program consists of a main routine and 13 subroutines. The 
main routine calls the subroutines in proper sequence and also provides 
for a sequential variation of the parameters r^, and V^. The subroutines 
and the function each performs are listed as follows: 



Subroutine INPUT2 


Reads in data cards 


Subroutine READY 2 


Calculates and tabulates energies, wave numbers, 
scattering lengths, effective ranges, potential 
parameters, and calls bessel functions. 


Subroutine FIGURE 


Prepares for and calls Integration routine, forms 
the logrithmic derivative of the Integrated function. 


Subroutine RKUTTA 


Integrates equation provided by DERIV for s 0. 


Subroutine RKUTTA 1 


Integrates equation provided by DERIVl for : 1. 


Subroutine DERIV 


Provides a 0 equation for Integration for a 
particular potential. 
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Subroutine 0ERIV1 

Subroutine FINISH 

Subroutine OUTPUT 

Subroutine BESSEL 
Subroutine SKETCH 

Subroutine APPEND 

Subroutine NOW 



Provides j( s I equation for integration for a 
particular potential. 

Computes phase shifts, cross sections, and 
normalized cross sections. 

Prints out resulting r^/a ratio, Vq, r^, energy, 
phase shifts, and cross sections. 

Calculates bessel functions 

Prepares for and calls DRAW to plot all cross 

sections on one graph for comparison. 

Stores base cross section for normalization 
in FINISH. 

Prepares for and calls DRAW to plot all norma- 
lized cross sections on one graph for comparison. 
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APPENDIX II 





PROGRAM DICTIONARY 


EINT 


Initial value of center of mass energy in MeV. 


EINC 


Incremental value of center of mass energy in MeV. 


1 1 


Maximum number of center of mass energy values. 


ROS 


Base effective range In cm. 


AS 


Base scattering length in cm. 


CHANGE 


Percentage incremental change in base effective range 
or base scattering length in decimal form. 


MANY 


Number of 4* or - incremental charges to be made. 


REND 


Radius at which integrated solution is matched to 
exact solution for phase shift determination in cm. 


NUMBER 


Number of integration steps. 


MODE 


Print out control for program monitoring 

■ 0 for result print out only 
z 1 for minor progress report 

■ 2 for print out of all calculations 


III 


Energy subscript for wave function print out. 


MOMENTA 


5 0 X- 0 calculations only 

si j(s 1 calculations only 

« 2 both {s 0 and 1 calculations. 


FK 


Wave number, k 


FK2 


2 

Wave number squared, k . 


FMANY 


MANY in floating point. 


RA 


Ratio of effective range to scattering length, r^j/a. 


S 


Well depth parameter, s 
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B 


Intrinsic range, b 


POT 


Well depth, Vq 


RN 


Well range, r^ 


ARGB 


Bessel function argument, kr 


VALUE 


Output array of bessel functions 


VAUO 


Spherical * 0 bessel function, j 


VAUl 


Spherical j|s 1 bessel function, J] 


VALJ2 


Spherical (s 2 bessel function, J2 


VALNO 


Spherical (>0 neuimian function, Hq 


VALNt 


Spherical ^ = 1 neumann function, n^ 


VALN2 


Spherical ^ 2 neumann function, n2 


A 


A particular potential range from array, RN 


V 


A particular potential depth from array, POT 


FNUMBER 


Number of integration steps in floating point 


STEP 


integration step size 


R 


Array of radius values at each integration step 


YVARS(l) 


Wave function 


WARS (2) 


Slope of wave functions 


XVAR 


Radius variable of integration 


Y 


Wave function at end value of integration step for j[ s 0 


YDOT 


Slope of wave function at end value of integration 
step for 11 ■ 0. 


COMPARO 


Logrithmic derivative of wave function at end of 
integration range for { s 0. 


COMPARl 


Logrithmic derivative of wave function at end of 
integration range for 1. 


XC 


Runge kutta weighted radius step 
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YC 


Dependent variable for XC 


DY 


Derivative of dependent variable for XC 


YL 


Same as Y for s 1 . 


YLOOT 


Same as YDOT for jf s 1 . 


DELTAO 


5 : 0 phase shift. 


DELTA I 


j| s 1 phase shift. 


XSECO 


^ s 0 cross section 


XSECI 


^ a 1 cross section 


TOTAL 


Sum of singlet and triplet cross sections 


BETTER 


Normalized ^=0 cross sections 


FNORM 


Base cross section to be used for normalization 


NUMBERR 


Number of simultaneous first order differential equa 




tions to be integrated. 
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APPENDIX III 
PROGRAM LISTING 



6h- 



MAIN PROGRAM 

DIMENSION E( 100) ,FK( 100) »FK2(100) ,RN( 100) .VALJO(IOO) .VALJK 100) 

1 ,VALJ2 ( 100) , VALNO(IOO), VALNl ( 100 ) »VALN2 ( 100 ) » R(IOOO), Y(IOOO), 
2YDOT(1000), COMPARO( 100) » COMPAR 1 ( 100 ) , YL ( 1 000 ) , YLDOT ( 1 000 ) , POT ( 1 
300) .DELTAO( 100) ♦DELTAl ( 100 ) ,XSEC0( 100) »XSEC1 (100) ♦TOTAL! 100) ♦RA ( 10 
40) ,FNORM(100) , BETTER ( 10 ♦ 100 ) 

COMMON E^EINT^EINCt FK^ FK2t RINT, RINCt RN, A^ II, LL , NN , MODE, 
1REND,VALJ0,VALJ1 ,VALJ2 ♦VALNO , VALNl ,VALN2 , JJ,M, NUMBER, R, I , L,V,K, 
2N,Y, YDOT, YL^ YLDOT, III, COMPARO^ COMP AR 1 , POT , V I NT , VI NC , MOMEN TA , 
3DELTA0^DELTA1 ,XSEC0,XSEC1 , TOTAL, RA,ROS^ AS ♦Change, MANY ♦$, 8, FNORM, BE 
4TTER 

CALL INPUT2 
CALL READY2 ■ 

A=RN(MANY+1) 

V=POT (MANY+1 )/. 41499 E-24 

CALL FIGURE 

CALL FINISH 

CALL APPEND 

DO 1 L=1,LL 

A=RN(L) 

V= POT( L) /.41499E-24 

CALL FIGURE , 

CALL FINISH 
CALL 0UTPUT2 

CALL SKETCH (E,XSECO, MOD, NUMPTS, I I ,L,LL) 

1 CONTINUE 

CALL NpW(E,BETTER,II,L,LL) 

END 
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SUBROUTINE INPUT2 

DIMENSION E( 100) ,FK( 100) »FK2(100) ^RNdOO) tVALJOdOO) tVALUldOO) 
l»VALJ2dOO), VALNOdOO). VALNl (100 ) ♦ VALN2 (100 ) * R(1000)» Y(1000)» 
2YDOT(1000), COMPARO( 100) . COMP ARl ( 100 ) » YL ( 1 000 ) tYLOOT (1000 ) »POT (1 
300) ♦DELTAOdOO) , DELT A1 ( 100 ) , XSECO ( 100 ) .XSECl ( 100) .TOTAL (100) »RA(10 
AO) ,FNORM (100) .BETTER (10.100) 

COMMON E.EINT.EINC. FK. FK2. RINT. RINC. RN, A. II, LL, NN. MODE.' 
1REND.VALJ0,VALJ1,VALJ2 .VALNO. VALNl .VALN2.JJ.M, NUMBER.' R. I .L,V.K. 
2;N.Y. YDOT, YL, YLDOT. III. COMPARO. COMPARl. POT, VINT .VINC. MOMENTA. 
3DELTA0,DELTA1 .XSECO, XSECl .TOTAL, RA.ROS, AS. CHANGE. MANY. S.B. FNORM .be 
ATTER 

READ l,EINT,EINC,ir 

1 FORMAT( 2E20.A. 110) 

READ 2, ROS. AS, CHANGE, MANY 

2 FORMAT( 2E20.6,F10,A, 110) [ 

READ 3, REND, NUMBER 

3 FORMAT (E20. A, 110) 

READ A, MODE. 

A ,F0RMAT( 110) • 

READ 5,III 
5 FORMAT ( 110) 

READ 6, MOMENTA 
6' FORMAT ( 110) 

IF(MODE - 1) 8,8,7 

7 PRINT 1,EINT,EINC,II 
PRINT 2, ROS, AS, CHANGE, MANY 
PRINT 3, REND, NUMBER] 

PRINT A, MODE 

PRINT 5, 1 1 I 
PRINT 6, MOMENTA 

8 RETURN 

end . . 
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SUBROUTINE FIGURE 

DIMENSION E( 100 ) ,FK( 100) »FK2 ( 1 00 ) , RN ( 100 ) .VALJOC 100) tVALJl (100) 

1 *VALJ2( 100) , VALNO(IOO), VALNl ( 100 ) » VALN2 UOO ) ♦ R(1000)» Y(1000)» 
2YDOT(lOOO), COMPARO( 100) ♦ COMPARl ( 100 ) » YL ( 1000 ), YLDOT ( 1000 ) »POT ( 1 
300) »DELTAO( 100) .DELTAl (lOO) ,XSECO (100) »XSEC1 (100) » TOTAL ( 100) »RA( 10 
40) ,FNORiM( 100) » BETTER ( 10»100) 

DIMENSION YVARSOO) 

COMMON E»EINT»EINC» FK» FK2» RINT» RINC» RN, A» II, LL, NN, MODE, 

irend,valjo,valji,valj2»valno»valni,valn2*JJ»m, number, R,I ,L,V,K, 

2N,Y, YDOT, YL» YLDOT, III» COMPARO, COMPAR 1 , POT , V I NT ♦ V INC ,MOMENTA » 
3DELTA0,DELTA1 ,XSEC0, XSEC1» TOTAL, RA,ROS, AS » CHANGE, MANY, S,B,FN0RM, BE 
4TTER 

.FNUMBER = NUMBER 
STEP = REND / FNUMBER 
R( 1 ) = 0. 

NN = NUMBER +1 r ■ 

DO 299 N=2,NN 
299 R(N) = R(N-l) + STEP 
DO 292 K=2,II 
■ YVARS(l) =0. 

YVARS(2) = 1. 

IF(MOMENTA - 1) 283,284,283 

283 CALL RKUTTA(2,XVAR,YVARS,STEP) 

IF(K-III) 293,297,293 

297 PRINT 296 

296 FORMAT (13X,7H ENERGY , 13X , 7H RADIUS, 6X, 14H WAVE FUNCTION) 

PRINT 295, E(K) 

295 FORMAT( F20.2 ) 

PRINT 294, (R(N), Y.(N), N = 1,NN) 

294 FORMAT ( E40 .4 ♦ E20 .4 ) 

293 COMPARO (K) =YD0T(NN-1 ) /Y(NN-1 ) 

284 IF (MOMENTA - 1) 2^2,291,291 

291 YVARS(l) = 0. 

YVARS(2) =1. 

CALL- RKUTTAl ( 2 ,XVAR , YVARS,STEP ) 

IF(K-III) 285,289,285 
289 PRINT 288 

288 FORMAT ( 13X , 7H ENERGY , 1 3X , 7H RADIUS, 6X, 18H WAVE FUNCTION L=l) ‘ 
PRINT 287, E(K) 

287 FORMAT( F20.2) 

PRINT 286, (R(N) ,YL(N) ,N = 1,NN) 

286 FORMAT( E40.4,E20.4) 

285 COMPARl (K) = YLD0T(NN-1 ) /YL(NN-1 ). 

292 CONTINUE 
RETURN' 

END 



64 , 



SUBROUTINE RKUTTA ( NUMBERR , X VAR , WARS » STEP ) 

DIMENSION E( 100 ) »FK( 100) »FK2 ( 100) ,RN( 100) .VALJ0( 100 ) tVALUK 100 ) 

1 ,VALJ2( 100) , VALN0(100)» VALNl ( 1 00 ) ♦ VALN2 ( 100 ) » R(IOOO), Y(1000)» 
2yDOT(1000), COMPARO( 100) » COMPARl ( 100 ) , YL ( 1000 ) , YLDOT ( lOOO ) »POT ( 1 
300) »DELTA0( 100) ,0ELTA1 ( 100 ) ,XSEC0( 100) *XSEC1 (100) ,T0TAL( 100) »RA( 10 
40) ,FNORM( 100) , BETTER (1 0 ♦ 100 ) 

DIMENSION YVARS(30)»AK(4,30)»DY(30) ,YC(30) *0(4) 

COMMON E*EINT,EINC» FK, FK2» RINT* RINC» RN , A» II* LL , NN , MODE* 
1REND*VALJ0*VALJ1 *VALJ2* VALNO* VALNl .VALN2 * JJ*M* NUMBER* R* I *L*V*K* 
2N*Y* YDOT* YL* YLDOT* III* COMPARO* COMPARl * POT *V I NT *V I NC * MOMENTA * 
3DELTA0*DELTA1 *XSECO*XSEC1*TOTAL*RA*ROS*AS*CHANGE*MANY*S*B*FNORM*BE 
4TTER 
C( 1 )=0.0 
C(2)=0.5 
C(3)=0.5 
C(4)=1.0 
DO 298 N=1*NN 
XVAR = R(N) 

DO 1 1=1*4 
XC=XVAR + C( I )*STEP 
bo 2 J=1*NUMBERR 

2 YC(J) =YVARS(J) + C( I )*AK(I-1*J) 

,CALL DERIV(XC*YC*DY) 

,DO 1 J=1*NUMBERR 
1 AK( I *J) = STEP * DY( J) 

DO 3 J = 1* NUMBERR 

3 YVARS(J) =YVARS(J) '+ ( A)C (1 * J ) +2 .*A)C ( 2 » J ) +2 .*AK ( 3 * J ) +AK ( 4* U ) ) /6 . 
;YD0T(N) = YVARS(2) 

298 Y(N) = YVARS(l) 

RETURN 

END 
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SUBROUTINE RKUTTAl (NUMBERR .X VAR .WARS ♦STEP ) 

DIMENSION E( 100) ,FK( 100) ♦FK2( 100) ,RN(100) tVALJOI 100) tVALJ! (100) 

1 ♦VALJ2( 100) ♦ VALN0(100)f VALNl ( 100 ) ♦ VALN2 ( 100 ) ♦ R(IOOO)^ Y(IOOO)^ 
2YDOT(1000)f COMPARO( 100) ♦ COMPARl ( 100 ) ♦ YL ( 1000 ) ♦YLDOT ( 1000 ) »POT (1 
300) ,DELTAO( 100) .DELTAl (100 ) ,XSECO( 100) ♦XSECl ( 100 ) ♦TOTAL ( 100 ) ♦RAdC 
40) ♦FNORM( 100) ♦BETTER ( 10^100 ) 

'DIMENSION YVARSOO) ,AK(4^30) ♦DY(30) ♦YC(30) ♦C(A) 

COMMON E^EINT^EINC^ FK^ FK2^ RINT^ RINC^ RN ♦ Ar !!♦ LL^ NN ♦ MODE^ 
i rend ♦VALJ 0^VALJ1,VALJ2^VALN0^ VALNl ♦VALN2^ JJ^M, NUMBER, R^ I ♦L^V^K^ 
2N,Y^ YDOT^ YL^ YLDOT ♦ Ild COMPARO ♦ COMP AR 1 ♦ POT ♦ V I NT ♦ V I NC ♦ MOMENTA , 
3DELTA0,DELTA1^XSEC0,XSEC1 ♦TOTAL ♦RA^ROS, AS, CHANGE, MANY, S^B^FNORM ♦BE 
4TTER 
C(1)=0.0 
C(2)=0«5 

C(3)=0.5 ■; 

C(4)=1.0 
DO 290 N=1,NN 
XVAR = R(N) 

DO 1 1=1,4 

:XC=XVAR + C( I )*STEP ' ' 

;DO 2 J=1,NUMBERR , 

2 YC(J) =YVARS(J) + C( I )*AK( 1-1, J) 

/CALL DERIVl (XC,YC,DY) 

iDO 1 J=1,NUMBERR 
1 AK( I ♦ J) = STEP * DY( J) . 

DO 3 J = 1, NUMBERR ' . 

3 YVARS(J) =YVARS(J) + ( AK ( 1 ♦ J ) +2 .*AK ( 2 ♦ J ) +2 .*AK ( 3 ♦ J ) +AK ( 4, J ) ) /6 . 
;YL(N) = YVARS(l) 

290iYLD0T(N) = YVARS(2) 

RETURN 

END 
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SUBROUTINE FINISH 

dimension E(100),FK( 10 0) »FK2(100) ,RN(100) »VALJ0(100) ♦VALJK 100) 

1 »VALJ2( 100) , VALNOdOO)* VALNl ( 100 ) ♦ VALN2 ( 100 ) ♦ R(1000)» Y(IOOO)* 
2YD0T(1000), COMPARO( 100) » COMPARl ( 100 ) ♦ Yl, { 1000 ) ,YLDOT ( 1000 ) tPOT ( 1 
300) tDELTAOl 100) »DELTA1 (100) ,XSECO( 100) ,XSEC1 (100) » TOTAL ( 100) .RA(10 
40) ,FNORM( 100) » BETT ER ( 1 0 » 100 ) 

COMMON E»EINT»EINC, FK, FK2» RINT, RINC» RN, A» II» LL, NN , MODE»’ 
1REND*VALJ0»VALJ1,VALJ2 »VALN0»VALN1 »VALN2 ♦ JJ*M, NUMBER* R» I ♦ L*V,K, 
2N*Y* YDOT* YL» YLDOT ♦ III* COMPARO* COMPAR 1 * POT * V I NT * V I NC * MOMENTA * 
3DELTA0*DELTA1 *XSEC0*XSEC1 * TOTAL* RA*ROS*AS* CHANGE *MANY*S*B*FNORM*BE 
4TTER 

IF (MOMENTA - 1) 399'* 397 * 399 
399 DO 398 1=2*11 

DELTAO( I ) =ATANF( ( REND* ( COMPARO ( I )*VALJO( I )+FK( I )*VALJ1 ( I ) )- 
1VALJ0( I ) ) /( REND* (COMPARO ( I )*VALNO( I )+FK( I )*VAlN1 ( I) ) - 
2VALN0( I ) ) ) 

398 XSECO( I )=4.*3.1416/FK2 ( I ) *S I NF ( DELT AO ( I ) )**2 
397 IF(M0MENTA-1 ) 394*396*396 
396 DO 395 I=2*II 

DELTAl ( I ) =ATANF( ( REND* (FK( I ) /3.*VALJ0( I )-2./3.*FK( I ) * - 

1VALJ2( I ) -COMPAR 1( I )*VALJ1( I ) )+VALJl( I ) ) /(REND*(FK( I ) /3.* 

2VALN0( I )-2./3.*FK( I ) *VALN2 ( I )-COMPARl(I )*VALN1(I) )+VALNl( I ) ) ) 

395 XSECK I )=12.*3.1416/FK2( I ) *S INF ( DELTAl ( I ) )**2 
390 DO 389 1=2*11 

389 TOTAL(I) = .25 * XSECO ( I ) + .75 * XSECKI) 

394 DO 900 1 = 2* 1 1 .f 

BETTER(L*I) = XSECO ( I ) /FNORM ( I ) 

90 0 CONTINUE 

RETURN 

END 
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SUBROUTINE OUTPUT2 

pi MENS I ON E( 100) »FK( 100) ♦ FK2 ( 100 ) » RN ( 100 )'* VAL JO ( 100 ) ♦VALJK 100) 
1»VALJ2(100) , VALN0(100)» VALNl ( 100 ) » VALN2 ( 100 ) » R(1000)» Y(IOOO). 
2YD0T(1000), C0MPAR0( 100)^ COMPARI ( 100 ) » YL( 1000 ) ,YLDOT ( 1000 ) »POT ( 1 
300) »DELTA0( 100) »DELT A1 ( 100 ) , XSECO ( 100 ) , XSECl (100) » TOTAL (100) ,RA(10 
40) »FNORM( 100) » BETTER (10,100) 

COMMON E,EINT,EINC, FK, FK2, RINT, RING, RN , A, II, LL, NN, MODE, 
1REND,VALJ0,VALJ1,VALJ2,VALN0,VALN1,VALN2, JJ,M, NUMBER, R, I ,L,V,K, 
2N,Y, YDOT, YL, YLOOT, III, COMPARO, COMPARI , POT ,V I NT , V I NC , MOMENTA , 
3DELTA0,DELTA1 , XSECO, XSECl, TOTAL, RA,ROS, AS, CHANGE, MANY, S,B,FNORM, BE 
4.TTER 

PRINT 600,RA(L) 

600 F0RMAT(34H EFFECTIVE RANGE/SCATTERING LENGTH ,/, E34 . 10 ) 

PRINT 599,A,POT(L) 

599 FORMAT (13X,7H RAD I US , 4X , 16H POTENTIAL DEPTH , / , E 20 .4 , E 20 . 4 ) 

PRINT 598 

598 F0RMAT(7H ENERGY, 12H PHASE ANGLE, 15H CROSS SECT I ON , 12H' PHASE ANGLE 
1,15H CROSS SECTION, 2X,6H TOTAL ,/, 1 IX ,4H L=0,7X,4H L=0,10X,4H L=l, 
2:8X,4H L = 1,8X,6H CROSS ,/, 63X , 8H SECTION) 

PRINT 597, ( E( I ) ,DELTA0 ( I ) , XSECO ( I ) ,DELTA1 ( I ) ,XSEC1 ( I ) , TOTAL ( I ) , 1=2 
1,11) 

597 F0RMAT(F7.2,E12.4,E15.4,E12.4,E15.4,E11.4) 

PRINT 596 • 

596; F0RMAT(7H ENERGY , 5X , 25H NORMALIZED CROSS SECTION) 

PRINT 595, ( E( I ) ,BETTER (L,I ) , 1=2, 1 1 ) 

595 FORMAT(F7.2,E20.8) • i 

RETURN 
END 
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SUBROUTINE BESSEL ( X . N . VALUE ♦ I NDEX , MM ) 

DIMENSION VALUE (40) ,BJ(100) ,F{ 100) 

IF (X) 10tl0,9 
9 IK=X/2.0 
K=N+20+IK 

IF(K-IOO) 300*299,299 

299 PRINT 298 

298 F,0RMAT(66H ARGUMENT PLUS ORDER TOO LARGE. INCREASE DIMENSIONS AND 
IRUN AGAIN.) 

RETURN 

300 CONTINUE 
Z=1.0/X 

PI=3. 141592653 

GOTO (40, 40, 60 *60, 80 ^90 ,100,110) , I NDEX 

10 PRINT 8,X 

80FORMAT ( IX , 63HSUBROUT I NE BESSEL DOES NOT WORK FOR X=0 OR LESS THAN 
1: 0, HERE X= E15.8) 

RETURN 

40 BJ(K+1)=0.0 , 

BJ(K )=10.0E-30 
L = K 

11 FL=L 

BJ( L-1 ) =2.0*FL*Z»BJ ( L)-BJ('l+1 ) 

IF(L-l) 12,12,13 
. 13 L = L-1 

GO TO 11 

12 SUM=0.0 
IK = K 

DO 14 1=2, IK, 2 

14 SUM=SUM+BJ( I ) 

C=1.0/(BJ(0 )+2.0*SUM) 

GO TO( 42,43 ) , INDEX , 

43 DO 45 I=0,K 
45 BJ( I )=C*BJ( I ) 

GO TO 54 
42 DO 41 1=0, N 
• 41 VALUE! I )=C»BJ( I ) 

RETURN 
50 GO TO 40 
54 EC=. 5772156649 
SUM=0.0 
IK=K/2 

DO 15 1 = 1, IK 
FI=I 

15 SUM=SUM+( (-1.0)**( I-1)*BJ(2*I ) ) /FI 

F(0)=( 2.0/PI )*( BJ(0)*( LOGF(X/2.0)+EC)+2.0*SUM) 

,F(1) = (BJ( l)*F(0)-2.0/(PI*X) )/BJ(0) 

IN=N-1 

DO 16 1=1, IN 
,FI = I 
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16 F( I+1)=2.0*FI*Z*F( I )-F( I-l ) 

DO 51 1=0, N 

51 VALUE( I )=F( I ) 

RETURN 

60 ' BJ(K+1 ) =0. 

BJ(K)=10.E-30 
' L=K 
19 FL=L 

B J ( L-1 ) =2 .*FL*Z^^BJ ( L ) +BJ { L + 1 ) 
IF(L-1)17,17,18 
18 L=L-1 
GO TO 19 

17 SUM=0.0 

DO 21 1=1, K 

21 SUM=SUM+BJ( I ) 

C=(EXPF(X) )/(BJ(0)+2.*SUM) 

GO T0( 1,1,61,62) , INDEX 

1 PRINT 1000 
RETURN 

lOOo' FORMAT(40H CARD MISSING. IGNORE RESULTS. 
62, .DO 64 1=0,1 

64 BJ( I )=C*BJ< I ) 

;.G0 TO 71 

61 DO 65 1=0, N 

65 VALUE! I )=C*BJ( I ) 

RETURN 

70 GO TO 60 

71 IF(5.0-X) 120,121,121 

121 M=10 

GO TO 122 
120 M=5 

122 F (0)=GAUSS(0.0,1.0,M,X) 
F(1)=(1./X-F(0)*BJ(1) )/BJ(0) 

IN = N-1 , -■ 

DO 22 1 = 1, IN 
F I = I 

22 F ( I+1)=2.0*FI*Z»F (I)+F (I-l) 

■ DO 72 1=0, N 

72. VALUE! I ) = F ! I ) 

^ RETURN 
80, ,F!K+1) =0.0 

F !K)=10.0E-30 
. L=K 
25 FL=L 

F !L-1)=!2.0*FL+1.0)*Z*F !L)-F !L+1) 

IF !L-1) 23,23,24 

- 24 L=L-1 

• GO TO 25 

23 C=SINF!X)/!X*F !0) ) 

IF!MM) 200,200,201 

201 ,SQ = SQRTF!2.*X/PI ) 

DO 202 1=0, N 

202 :VALUE! I )=C*SQ*F !I) 

RETURN 
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200 DO 26 I=0,N 

26 VALUE( I )=C*F ( I ) 

RETURN 

90 F (0)=-COSF(X) /X 

F ( 1)=-SINF(X) /X-COSF(X)/X**2 
IN=N-1 

DO 27 1 = 1, IN 
FI = I 

27 F( I+1)=(2,*FI+1.0)*Z*F(I )-F( I-l) 

IF(MM) 203,203,204 

204 SQ=SQRTF(2.*X/PI ) 

DO 205 1=0, N 

205 VALUE(I)=F ( I ) *SQ* ( -1 . ) ** ( - I "1 ) 

RETURN 

• 203 DO 91 1=0, N 

91 VALUE( I )=F ( I ) 

* RETURN 

100 F (K+1)=0.0 

F (K)=10.0E-30 . , 

30 FL=L 

F ( L-1 ) =( 2.0*FL+1,0) *Z*F ( L ) +F (L+1) 

IF (L-1) 28,28,29 • 

29.L=L-1 
GO TO 30 

28 C=SINH (X ) / (X*F ( 0 ) ) 

IF(MM) 206,206,207 

207 SQ=SQRTF( 2.*X/PI ) 

DO 208 1=0, N 

208 VALUE( I ) =C*SQ*F (I) 

RETURN 

206 DO 101 1=0, N 

101 VALUE( I )=C*F ( I ) ■ ■ 

RETURN 

110 F (0)=-PI*EXPF(-X)/(2.0*X) 

f (1)=(PI*(X+1.0)*EXPF(-X) )/(2.0*X**2) 

•IN = N-1 

DO 31 1 = 1, IN 
FI = I 

31 F (I+1)=F ( I-1)-(2.0*FI+1.0)*Z*F (I) 

IF(MM) 209,209,210 

210 SQ=SQRTF(2.*X/PI ) 

DO 211 1=0, N 

211 VALUE(.I)=F ( I )*SQ*(-1. )*^^(-I-l ) 

RETURN 

209 DO 111 1=0, N 

111 ;VALUE( I ) = F ( I ) 

■RETURN 
'END 

FUNCTION GAUSS NUMERICAL INTERGRATION 
, FUNCTION GAUSS REQUIRES FUNCTION GAUSSF(X) NUMERICAL INTERGRATION 
FUNCTION GAUSS(G,B,M,X ) 

•DIMENSION R( 10) ,U( 10) 

,U( 1 )=+. 0744371695 
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iU(2)=-.074A371695 
U( 3 )=+. 2166976971 
U(4)=-. 2166976971 
U(5)=+. 3397047841 
U(6)=-. 3397047841 
U(7)»+, 4325316833 
U( 8 )*-. 4325316833 
U(9)=+. 4869532643 
U( 10)=-*4869532643 
R( 1 )=+. 1477621124 
R(2)=+. 1477621124 
:R( 3 )=+. 1346333597 .. 

R(4)=+. 1346333597 
R(5)=. 10954318113 
R(6)=. 10954318113 
R(7)=+. 07472567458 
R(8)=+. 07472567458 
R(9)=+. 03333567215 
R( 10)=+. 03333567215 
A=G 
FN=M 

P=(B-A)/FN 
GAUSS=0.0 
DO 1 J=1*M 
C=0.0 
Y=A+P 

DO 2 1=1.10 

D=(R( I )*GAUSSF( ( Y-A)*U( I ) + ( A+Y ) /2 . 0 , X ) ) * ( Y-A ) 
2 C=D+C 
A = Y 

1 GAUSS=GAUSS +C 
RETURN 
END 

FUNCTION GAUSSF(U.X) 

IF ( 1.0 /U-710.0) 123.124,124 
124 GAUSSF=0.0 
RETURN 

123 ‘gAUSSF= (EXPF(-X*COSH(1.0/U-1.0) ) )/U**2 
END 
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SUBROUTINE SKETCH ( E » XSECO » MOD .NUMPTS » I I *L,LL) 

DIMENSION E( 100) , XSECO { 100 ) , XSECl (100) .TOTAL (100) 

DIMENSION X(900) ,Y(900) ,ITITLE( 12) 

LABEL = 4H 
DO 1 1=1.12 

1 ITITLE( I ) = 8H 
ITITLE( 1 ) =8HOSGOOD.F 
ITITLE( 2) =8H.B. .BOX 
ITITLE(3)=8H09 
ITITLE(4) =8HEXPONENT 
ITITLE( 5)=8HIAL WELL 
ITITLE(6)=8H RUN 15 
ITITLE(7)=8H GRAPH 0 
ITITLE(8)=8HF CROSS 
ITITLE(9)=8HSECTION 
ITITLE( 10) =8HVS ENERG 
ITITLE( 11 )=8HY FOR L= 

ITITLE(12) = 8H0 
NUMPTS = II-l ; 

IF(L-1)2.2.3 

2 MOD = 1 
GO TO 6 

3 IF(LL-L) 5.5.4 

4 MOD = 2 
GO TO 6 

5 MOD = 3 

6 CONTINUE 

CALL DRAW ( NUMPTS. E.XSECO. MOD. 0. LABEL. I TI TLE. .3.0. O.O.O.O. 9 . 15.1 . 
lAST) 

RETURN 

END 
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SUBROUTINE APPEND 

DIMENSION E ( 100) »FK( 100) .FK2 ( 100 ) , RN ( 100 ) iVALJ0{ 100) »VALJ1 (100) 
1»VALJ2(100), VALNO(IOO), VALNl ( 100 ) »VALN2 ( 100 ) ♦ R(1000)» Y(1000)» 
2YDOT(1000), COMPAROdOO)* COMPAR 1 ( 1 00 ) » YL ( 1 000 ) , YLDOT ( 1 000 ) tPOT ( 1 
300) »DELTAO( 100) ,DELTA1 (100 ) ,XSECO( 100) ,XSEC1 ( 100) »TOTAL{ 100) »RA (10 
40) »FNORM( 100) ♦ BETTER ( 10 ♦ 100) 

COMMON E»EINT»EINC» FK» FK2» RINT» RING, RN, A, II, LL, NN, MODE, 
1REND,VALJ0,VALJ1 ,VAUJ2 ,VALNO, VALNl ,VALN2 ,JJ,M, NUMBER, R, I , L,V,K, 
2N,Y, YDOT, YL, YLDOT, III, COMPARO, COMPAR 1 , POT ,V I NT ,V I NC , MOMENTA , 
3DELTAO,DELTA1,XSECO,XSEC1,TOTAL,RA,ROS,AS,CHANGE,MANY,S,B,FNORM,BE 
4TTER 

DO I 1=2,11 
I FNORM( I ) = XSECO( I ) 

RETURN 

END 
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SUBROUTINE NOW ( E ♦ BET TER , I I » L , LL ) 

DIMENSION E(IOO), BETTER ( 1 0 » 100 ) , X ( 900 ) ♦ Y ( 900 ) ♦ I T I TLE (12 ) 

LABEL = AH 
DO 1 1=1 *12 
1 ITITLE(I)=8H 

ITITLE( 1)=8H0SG00D,F 
ITITLE( 2) =8H.B. ,BOX 
ITITLE( 3) =8H09 
ITITLE(7)=8HNORMALIZ 
ITITLE(8)=8HED CROSS 
ITITLE(9)=8H SECTION 
ITITLE( 10)=8H FOR L=0 
NUMPTS = II-l 
DO 7 L=1,LL 
IF(L-l) 2*2,3 
2. M0D=1 
GO TO 6 

3 IF(LL-L) 5, 5, A ' , 

A MOD = 2 ' 

GO TO 6 

5 MOD = 3 

6 CONTINUE 

DO 20 1=2,11 
20 BET( I )=BETTER(L,I ) 

CALL DRAW(NUMPTS,E,BET, MOD , 0 , LABEL', I T I TLE , . 3 , 0 , 0 , 0 , 0 , 0 , 9 , 15 , 1 , L 

lAST) 

7 CONTINUE 

RETURN ^ 

END 
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THE FOLLOWING SUBROUTINES ARE DEPENDENT UPON A PARTICULAR POTENTIAL SHAPE. 
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GAUSSIAN POTENTIAL 



SUBROUTINE READY2 

DIMENSION E(IOO) ,FK( 100) »FK2(100) »RN(100) »VALJ0(100) »VALJ1 (100) 

1 ,VALJ2( 100) ♦ VALN0(100)» VALNl ( 100 ) ♦ VALN2 ( 100 ) ♦ R(IOOO)* Y( 1000)» 
2YDOT(1000)» COMPARO( 100 ) ♦ COMPAR 1 ( 100 ) . YL ( 1 000 ), YLDOT ( 1000 ), POT { 1 
3 00) ,DELTA0( 100) .DELTAl (100 ) »XSEC0( 100) ,XSEC1 (100) » TOTAL ( 100) »RA( 10 
AO) ,FNORM(10 0) , BETTER ( 1 0 , 100 ) 

DIMENSION VALUE(AO) 

COMMON E,EINT»EINC, FK, F)C2» RINT* RING* RN, A* II, LL , NN , MODE, 

1 REND, VAL JO, VALJ 1 ,VALJ2,VALN0, VALNl ,VALN2 »JJ»M, NUMBER, R, I ,L,V,K, 
2N,Y, YDOT, YL, YLDOT, III, COMPARO, COMPARl , POT , V I NT , V I NC , MOMENTA , 
3DELTA0,DELTA1 ,XSEC0,XS EC 1 , TOTAL, RA,ROS, AS, change, MANY, S,B,FNORM, BE 
4TTER 

E( 1 ) = EINT 

DO 199 I = 2,11 

E( I ) = E( I-l ) + EINC 

FK(I) = SQRTF(E ( I ) /.41499E-2A) 

199 FK2 ( I )= FK( I ) * FK( I ) 

LL=2*MANY +1 
JJ = LL 

FMANY = MANY 

GAUSSIAN DISTRIBUTION 

RA(1) = ROS/AS*(l.- FMANY*CHANGE) 

DO 170 L=1,LL 

S=(1.“2.03*RA(L) + .58 7»RA(L)**2) / (1.-2.81*RA(L )+1.805*RA(L)*<^2) 
SS=S-1. 

B=RA(L)*AS/( l.-.641*SS+.568*SS**2-.721*SS**3+1.172*SS**4-1.015*SS* ' 
1*5 ) 

POT(L)=S*.229208E-23/B**2 

RN(L)=SQRTF(B**2/2.0604) 

170 RA(L+1) = RA(L) + CHANGE *ROS/AS 
IF(MODE “ 1) 190,192,197 

197 PRINT 196 

196 FORMAT (13X,7H ENERGY , 5X , 1 2H WAVE NUMBER , 3X , 20H WAVE NUMBER SQUARE 
ID) 

• PRINT 195 ,(E(I), FK( I ) ,FK2( I ) ,1=2, I I ) 

195 FORMAT (3E20.4) 

PRINT 171 

171 FORMAT(34H EFFECTIVE RANGE/SCATTERING LENGTH) 

PRINT 172,(RA(L) L=1,LL) 

172 F;ORMAT(E34.10) 

FORINT 194 

194 FORMAT (13X,7H RADIUS) 

PRINT 193, (RN(L) , L=1 ,LL) 

193 FORMAT (E20.4) 

PRINT 1192 

1192 FORMAT (16H POTENTIAL DEPTH) 

PRINT 2192, (POT(J) ,J=1,JJ) 

2192 FORMAT(F16.2 ) 

192 PRINT 191 
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191 FORMAK 34H RADIUS AND ENERGY VALUES COMPUTED) 
190 DO 187 I=2»II 

ARGB = FK( I ) * REND 

CALL BESSEL( ARGB ♦ 2 ♦ VALUE » 5 , 0) 

VALJOC I ) = VALUE! 0 ) 

VALJl ( I ) - VALUE! 1 ) 

VALJ2 ! I) = VALUE! 2 ) 

CALL BESSEL ! ARGB , 2 » VALUE» 6»0) 

VALNO! I ) = VALUE! 0 ) 

VALNl! I ) = VALUE! 1 ) 

VALN2! I ) = VALUE! 2 ) 

IF!MODE - 1) 187,187»189 
189 PRINT 188,1 
188 F0RMATI14I5) 

187 CONTINUE 

IF,!MODE - 1) 184,186,186 
186 PRINT 185 

185 FORMAT! 26H BESSEL FUNCTIONS COMPUTED) 

184 RETURN 
END 
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EXPONENTIAL POTENTIAL 



SUBROUTINE READY2 

DIMENSION E( 100) ,FK( 100) »FK2 ( 100 ) » RN ( 100 ) »VALJ0( 100 ) »VALJ1 ( 100 ) 

1 ,VALJ2 ( 100) , VALNOdOO)* VALNl ( 100 ) , VALN2 ( 100 ) » R(IOOO), Y(lOOO), 
2YD0T( 1000), COMPaRO( 100) , COMPARl ( 100 ) , YL ( 1 000 ) , YLDOT ( 1 000 ) » POT ( 1 
3 00 ) ,DELTAO( 100) »DELTA1 (100 ) ,XSECO (100) ♦XSECl (100) ♦ TOTAL ( 100 ) »RA ( 10 
40 ) ,FNORM( 100 ) , BETTER ( 10»100 ) 

DIMENSION VALUE(40) 

COMMON E,EINT»EINC, FK» FK2* RINT»’RINC, RN , A, II, LL, NN, MODE, 

1 REND »VALJ0,VALJ1»VALJ2»VALN0» VALNl »VALN2»JJ»M, NUMBER, R»I »L,V,K, 
2N,Y, YDOT, YL, YLDOT, III* COMPARO, COMPAR 1 , POT , V I NT , V I NC , MOMENTA , 
3DELTA0,DELTA1 ,XSECO, XSECl, TOTAL, RA,ROS, AS* CHANGE, MANY *5,6, FNORM* BE 
4TTER 

E(l) = EINT 

DO 199 I = 2,11 

Ed) = E(I-l) + EINC 

FKd) = SQRTF(Ed ) /.41499E-24) 

199 FK2( I )= FK( I ) * FK( I ) 

LL=2*MANY + 1 
JJ = LL 

EXPONENTIAL WELL 
FMANY = MANY 

RA( 1 )=ROS/AS*( 1.+FMANY*CHANGE) 

DO 170 L=1,LL 

S=(1.-2.9739*RA(L)+1.360 2*RA(L)^^*2)/ (l.-3.7134*RA(L)+3.1177*RA( D* 
1*2 ) 

SS=S-1. 

B=RA(L)*AS/ (1.-.904*SS+.745*SS**2-.542*SS**3+.174*SS**4 ) 

POT ( L) =S*.751541E-23/B**2 
RN(L)=B/3.5412 

170 RA(L+1) = RA(L) - CHANGE *ROS/AS 
IF(M0DE - 1) 190,192*197 

197 PRINT 196 

196 FORMAT (13X*7H ENERGY * 5X * 12H WAVE NUMBER * 3X , 20H WAVE NUMBER SQUARE 
ID) 

• PRINT 195 *(E(I)* FK( I ) *FK2 ( I ) * 1=2 * I I ) 

195 FORMAT (3E20.4) 

PRINT 171 

171 F0RMAT(34H EFFECTIVE RANGE/SCATTERING LENGTH) 

PRINT 172*(RA(L) L=1*LL) 

172 FORMAT (E34. 10) 

PRINT 194. 

194 FORMAT (13X*7H RADIUS) 

PRINT 193* (RN(D* L = 1.*LL) 

193 FORMAT (E20,4) 

PRINT 1192 

192 FORMAT (16H POTENTIAL DEPTH) 

PRINT 2192* (POT( J) *J=1*JJ) 

192 FORMAT(F16.2 ) 

192 PRINT 191 
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^ 191 FORMAK 34H RADIUS AND ENERGY VALUES COMPUTED) 

190 DO 187 1=2,11 

ARGB = FK( I ) * REND 

CALL BESSEL( ARGB , 2 , VALUE , 5 » 0) 

VALJO( I ) = VALUE( 0 ) 

VALJK I ) = VALUE( 1 ) 

VALJ2( I ) = VALUE( 2 ) 

CALL BESSEL ( ARGB , 2 , VALUE » 6*0) 

VALNO( I ) = VALUE( 0 ) 

VALNl ( I ) = VALUE( 1 ) 

VALN2 ( I ) = VALUE( 2 ) 

IF(MODE - 1) 187,187,189 
189 Pf^INT 188,1 
188 FORMAK 1415) 

187 CONTINUE 

IF(MODE - 1) 184,186,186 
186 PRINT 185 

185 FORMAT( 26H BESSEL FUNCTIONS COMPUTED) 

184 RETURN 
END 
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YUKAWA POTENTIAL 



SUBROUTINE READY2 

DIMENSION E{100),FK(100) »FK2 ( 100 ) , RN ( 1 00 ) ,VALJ0 { 100 ) »VALJ1 { 100) 

1 >VALJ2 ( 100) ♦ VALNO(IOO), VALNl ( 100 ) ♦ VALN2 ( 100 ) ♦ R(IOOO), Y(1000)> 
2YOOT(1000), COMPARO< 100) ♦ COMPARI { 100 ) ♦ YL ( 1000 )♦ YLDOT ( 1000 )♦ POT ( 1 
300) »DELTAO{ 100) »DELTA1 ( 100 ) , XSECO ( 100 ) ,XSEC1 ( 100) ♦ TOTAL (100) »RA(10 
40) »FNORM( 100) .BETTER ( 10.100) 

DIMENSION VALUE(40) 

COMMON E.EINT.EINC. FK . FK2. RINT. RING. RN. A. II. LL. NN. MODE. 

1 REND. VAL JO. VALJ1.VALJ2.VALN0. VALNl . VALN2.. J J . M . NUMBER. R. I .L.V.K. 
'2N.Y. YDOT. YL. YLDOT. III. COMPARO . COMPARI . POT .V I NT .V I NC . MOMENTA . 
3DELTAO.DELTA1.XSECO.XSEC1.TOTAL.RA.ROS.AS.CHANGE.MANY.S. 8 . FNORM.be 
4TTER 

E(l) = EINT 

DO 199 I = 2. 1 1 

E( I ) = E( I-l ) + EINC 

FK(I) = S0RTF(E( I )/.41499E-24) 

199 FK2 < I ) = FK< I ) * FK( I ) 

LL=2*MANY + 1 
JJ = LL 

FMANY = MANY 

RA( 1 )=ROS/AS*( l.+FMANY*CHANGE) 

C YUKAWA WELL 

DO 170 L=1.LL 

S= ( 1 .-3.5093-»RA(L)+2.0 354*RA(L) **2 ) / ( 1 . -“4. 1348*RA ( L ) +3. 7072 *RA ( L ) * 
1*2 ) 

SS=S-1. 

B=RA < L ) *AS/ { 1 «-l. 369*SS+1 . 093*SS**2-1 . 127*SS**3+1 .00 5*SS**4-.324*S 
1S**5 ) 

POT (L) =S*. 147585E-23*2.1196/B**2 
RN(L)=B/2.1196 

170 RA(L+1) = RA(L) - CHANGE *ROS/AS 

IF(MODE - 1) 190.192.197 

197 PRINT 196 

196 FORMAT (13X.7H ENERGY . 5X . 12H WAVE NUMBER . 3X . 2 0H WAVE NUMBER SQUARE 

ID) 

PRINT 195 .(Ed). FK( I ) .FK2( I ) . 1=2. I I ) 

195 FORMAT (3E20.4) 

PRINT 171 

171 FORMAT! 34H EFFECTIVE RANGE/SCATTERING LENGTH) 

PRINT 172. (RAIL) L=1.LL) 

172 FORMAT(E34.10) 

PRINT 194 

194 FORMAT (13X.7H RADIUS) 

PRINT 193. (RN(L) . L=1.LL) 

193 FORMAT (E20.4) 

PRINT 1192 

1192 FORMAT ( 16H POTENTIAL DEPTH) 

PRINT 2192. (POT( J) .J=1.JJ) 

2192 FORMAT! F16. 2 ) 
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N 192 PRINT 191 

' 191 FORMAT! 34H RADIUS AND ENERGY VALUES COMPUTED) 

190 DO 187 1=2,11 

ARGB = FK( I ) * REND 

CALL BESSEL! ARGB , 2 , VALUE , 5 » 0) 

VALJO! I ) = VALUE! 0 ) 

VALJl! I ) = VALUE! 1 ) 

VALJ2! I ) = VALUE! 2 ) 

CALL BESSEL ! ARGB ♦ 2 » VALUE , 6,0) 

VALNO! I ) = VALUE!0) 

VALNl ! I ) = VALUE! 1 ) 

VALN2 ! I ) = VALUE! 2 ) 

IF!MODE - 1) 187,187,189 
189 PRINT 188,1 . 

188 FORMAT!14I5) 

187 CONTINUE 

IF!MODE - 1) 184,186,186 
186 PRINT 185 

185 FORMAT! 26H BESSEL FUNCTIONS COMPUTED) 

184 RETURN 
END 
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SQUARE WELL POTENTIAL 



SUBROUTINE READY2 

DIMENSION E ( 100 ) » FK ( 100) ♦FKZ ( 100 ) » RN ( 100 ) » VAL JO ( 100 ) » VALJl ( 100 ) 

1 ,VALJ2( 100) » VALNO(IOO), VALNl ( 100 ) , VALN2 ( 100 ) » R(1000)» Y(1000)» 
2YDOT( 1000)» COMPARO ( 100 ) » COMP AR 1 ( 1 00 ) » YL ( 1 000 ) » YLDOT ( 1 000 ) ♦ POT ( 1 
300 ) »DELTAO (100) »DELTA1(100)»XSEC0(100) »XSEC1(100) »TOTAL(100) »RA{10 
40) ,FNORM( 100) , BETTER ( 10»100) 

DIMENSION VALUE(40) 

COMMON E»EINT»EINC» FK» FK2» RINT* RINC» RN, A» II* LL , NN, MODE* 
1REND*VALJ0* VALJl *VALJ2 *VALNO * VALNl * VALN2 * JJ»M* NUMBER* R* I * L*V* K* 
2N,Y* YDOT* YL* YLDOT* III* COMPARO* COMPAR 1 * POT * V I NT * V I NC * MOMEN TA * 
3DELT AO *DELTA1 »XSECO* XSECl * TOTAL *RA * ROS * AS * CHANGE * MAN Y * S * B * FNORM * BE 
4TTER 

E(l) = HINT 

DO 199 I = 2*11 

E( I ) = E( I-l ) + EINC 

FK(I) = SQRTF(E( I ) /.41499E-24) 

199 FK2( I )= FK( I ) * FK( I ) 

LL=2*MANY + 1 
JJ = LL 

SQUARE WELL 
FMANY = MANY 

CHANGE SIGN TO - FOR NORMAL OPERATION 
RA( 1 )=ROS/AS*( l.+FMANY*CHANGE) 

D0.170 L=1*LL 

S=( l.-1.574*RA(L)+.353*RA(L)**2)/( 1.-2.384*RA(L)+1.30*RA(L)**2) 

SS=S-1* 

B=RA (L)*AS/ ( 1 .-,5*SS+.486*SS**2-.476*SS**3+.378*SS**4-.158*SS**5) 
POT(L) = S*.102276E-23/B**2 
RN(l-) =B 

170 RA(L+1) = RA(L) - CHANGE *ROS/AS 
IF(MODE - 1) 190*192*197 

197 PRINT 196 

196 FORMAT (13X*7H ENE;RGY * 5X * 1 2H WAVE NUMBER * 3X,* 20H WAVE NUMBER SQUARE 
ID) ■ ' . . 

PRINT 195 *(E(D* FK ( I ) * FK2 ( I ) * I =2 , 1 1 ) 

195' FORMAT (3E20.4) 

PRINT 171 

171 FORMAT(34H EFFECTIVE RANGE/SCATTERING LENGTH) 

PRINT 172*(RA(L) L=1*LL) 

172 FORMAT! E34. 10) 

PRINT 194 

194 FORMAT (13X*7H RADIUS) 

PRINT 193* (RN(D* L=,1»LL) 

193 FORMAT (E20.4) 

PRINT 1192 ; 

1192 FORMAT ( 16H POTENTIAL DEPTH) 

PRI.NT 2192* ( POT( J ) * J = 1 *JJ) 

2192 FORMAT(F16.2 ) 

192 PRINT 191 
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191 FORMAK 34H RADIUS AND ENERGY VALUES COMPUTED) 
190 DO 187 1=2,11 

ARGB = FK( I ) * REND 

CALL BESSEL( ARGB , 2 , VALUE , 5 , 0) 

VALJO( I ) « VALUE! 0 ) 

VALJK I ) = VALUE! 1 ) 

VALJ2 ! I ) = VALUE! 2 ) 

CALL BESSEL ! ARGB , 2 , VALUE , 6,0) 

' • VALNO! I ) = VALUE!0) 

VALNl! I ) = VALUE! 1 ) 

VALN2! I ) = VALUE!2 ) 

IF!MODE - 1) 187,187,189 
189 PRINT 188,1 
188 FORMAT! 1415) 

187 CONTINUE 

IF!MODE - 1) 184,186,186 
186 PRINT 185 

185 FORMAT! 26H BESSEL FUNCTIONS COMPUTED) 

184 RETURN 



END 
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GAUSSIAN POTENTIAL 



SUBROUTINE DERIV (XC,YC,DY) 

DIMENSION £(100)>FK(100)*FK2(100)»RN(100) »VALJO( 100 ) »VALJ1 ( 100 ) 

1 ,VALJ2(100) , VALNO(IOO)* VALNl ( 100 ) , VALN2 ( 100 ) » R(1000)» Y(1000)» 
2YDOT(1000), COiMPARO( 100) ♦ COMPARl ( 100 ) ♦ YL ( 1000 ), YLDOT ( 1000 ) *POT { 1 
300 ) ,DELTA0( 100) ,DELTA1 (100 ) ,XSEC0 ( 100) »XSECl (100) » TOTAL ( 100 ) >RA ( 10 
40) ,FNORM(100) . BETT ER ( 1 0 » 1 00 ) 

DIMENSION YC(30), DY(30) 

COMMON E»EINT*EINC» FK, FK2» RINT» RINC» RN» A» II, LL, NN, MODE, 

1 REND, VALJ0,VALJ1,VALJ2,VALN0, VALNl ,VALN2, JJ,M, NUMBER, R, l', L,V,K, 
2N,Y, YDOT, YL, YLDOT , III, COMPARO, COMPARl , POT ,V I NT ,V I NC , MOMENTA , 
3DELTA0,DELTA1 ,XSEC0,XSEC1, TOTAL, RA,ROS,AS,CHANGE, MANY, S,B,FNORM, BE 
4TTER 

DY(1) =YC(2) 

DY(2) = - (FK2(K) + V * EXPF(-(XC/A) ** 2)) * YC ( 1 ) 

RETURN 

END 



SUBROUTINE DER I VI ( XC , YC , DY ) 

DIMENSION E( 100) ,FK( 100) ,FK2(100) , RN( 100 ) , VALJO ( 100 ) ,VALJ1 (100) 

1 ,VALJ2( 100) , VALNO(IOO), VALNl ( 100 ), VALN2 ( 100 ) , R(IOOO), Y(IOOO), 
2YDOT(1000), COMPARO( 100) , COMPARl ( 100 ) , YL ( 1 000 ), YLDOT ( 1 000 ), POT ( 1 
300 ) ,DELTAO( 100) ,DELTA1 (100) ,XSECO( 100) ,XSEC1 (100) ♦ TOTAL ( 100) ,RA (10 
40) ,FNORM( 100) , BETTER ( 1 0 , 100 ) 

DIMENSION YC(30), DY(30) 

COMMON E,EINT,EINC, FK, FK2, RINT, RINC, RN, A, II, LL, NN , MODE, 

1 REND, VALJO, VALJ1,VALJ2 ,VALNO, VALNl ,VALN2, JJ,M, NUMBER, R, I , L,V,K, 
2N,Y, YDOT, YL, YLDOT, III, COMPARO, COMPAR 1 , POT , V I NT , V I NC , MOMENTA , 
3 DELT AO, DELTA 1, XSECO, XS EC 1 , TOTAL, RA,ROS, AS, change, MANY ,S,B,FNORM, BE 
4TTER 

DY( 1 ) = YC(2) 

DY(2) = -(FK2(K) + V * EXPF(-(XC/A)**2) - 2. /XC**2 ) *YC ( 1 ) 

RETURN 

END 



1 
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EXPONENTIAL POTENTIAL 



SUBROUTINE DERIV (XC,YC,DY) 

DIMENSION E( 100) ,FK( 100) »FK2 ( 1 00 ) » RN ( 100 ) »VALJO( 100) »VALJ1 (100) 

1 ♦VALJ2( 100) ♦ VALN0(100)» VALNl ( 100 ) * VALN2 ( 100 ) » R(1000)» Y(1000)» 
2YDOT(lOOO), COMPARO( 100) » COMPARl ( 100 ) ♦ YL ( 1 000 )♦ Y LOOT ( 1 000 ) »POT ( 1 
300 ) ,DELTA0( 100) ,DELT A1 ( 100 ) ,XSEC0(100) ,XSEC1 (100) ♦TOTAL ( 100 ) tRA ( 10 
40) ♦FNORM( 100) ♦ BETTER ( 10 ♦100) 

DIMENSION YC(30)^ DY(30) 

COMMON E^EINT^EINC^ FK* FK2, RINT^ RINC^ RN^ A^ !!♦ LL ♦ NN^ MODE^ 
1REND^VALJO^VALJ1^VALJ2^VALNO^ VALNl ♦VALN2^JJ^M^ NUMBER ♦ R^I ♦L^V^K^ 
2N^Y^ YDOT^ YL^ YLDOT ♦ HU COMPARO^ COMPARl ♦ POT ♦V INT ♦V I NC ♦ MOMENTA ♦ 
3DELTA0 ♦ DELTAl ♦XSECO ♦ XSECl ♦ TOTAL ♦RA^ROS ♦AS ♦CHANGE ♦MANY ♦S^B^FNORM^ BE 
4TTER 

DY ( 1 ) = YC( 2 ) 

DY(2) = - (FK2(K) + V * EXPF(- XC/A)) * YC(1) 

RETURN 

END 



SUBROUTINE DER I VI ( XC ♦ YC ♦ DY ) 

DIMENSION E ( 100) ♦FK( 100) ♦FK2 ( 100 ) ^RN( 100 ) ♦VALJO( 100) ♦VALJl (100) 

1 ♦VALJ2 ( 100) ♦ VALNO(IOO)^ VALNl ( 100 ) ♦VALN2 ( 100 ) ♦ R(IOOO)^ Y(IOOO)^ 
2YD0T(1000)^ C0MPAR0( 100) ♦ COMPARl ( 100) ♦ YL ( 1 000 )♦ Y LOOT ( 1 000 ) ♦POT ( 1 
300 ) ♦DELTAO( 1 00 ) ♦DE LT A1 ( 1 OO ) ♦XSEC0( 100) ♦XSECl (100) ♦ TOTAL ( 100 ) ♦RA ( 10 
40) ♦FNORM( 100) ♦BETTER ( lO^lOO) 

DIMENSION YC(30)^ DY(30) 

COMMON E^EINT^EINC^ FK^ FK2^ RINT^ RINC^ RN^ A^ II, LL^ NN ♦ MODE^ 
IREND^VALJO ♦VALJl ♦VALJ2^VALN0 ♦VALNl ♦VALN2S JJ^M^ NUMBER ♦ R^ I ♦L^V^K^ 
2N^Y^ YDOT^ YL^ YLDOT ♦ IIH COMPARO^ COMPARl ♦ POT ♦V I NT ♦VI NC ♦ MOMENTA ♦ 
3DELTA0^ DELTA l^XSECO^XSECl^ TOTAL ♦RA^ROS ♦ AS ♦ CHANGE ♦ MANY ♦ S ♦ B ♦ FNORM ♦ BE 
4TTER 

DY( 1 ) = YC( 2 ) 

DY(2) = - (FK2(K) + V * EXPF(- XC/A) -2. /XC**2 ) *YC ( 1 ) 

RETURN 

END 
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YUKAWA POTENTIAL 



SUBROUTINE DERIV (XC,YC»DY) 

DIMENSION E ( 100 ) ,FK( 100) »FK2 ( 100) ,RN( 100) ,VALJ0( 100) ,VALJ1 ( 100) 

1 *VALJ2 ( 100) , VALN0(100)» VALNl ( 100 ) , VALN2 ( 100 ) > R(IOOO), Y(1000)» 
2YDOT(1000), COMPAROI 100) » COMPARl ( 100 ) , YL ( 1 000 ) , YLDOT ( 1 000 ) , POT ( 1 
300) ♦DELTAO(IOO) .DELTAKIOO) ,XSEC0(100) ,XSEC1(100) ,TOTAL(100) ,RA(10 
40) ♦F.NORMI 100) , BETTER ( 1 0 » 1 00 ) 

DIMENSION YC(30), DY(30) 

COMMON E»EINT»EINC, FK, FK2 » RINT, RING, RN, A» II, LL, NN, MODE, 

1 REND, VAL JO, VALJl ,VALJ2 ,VALNO, VALNl ,VALN2 , JJ,M, NUMBER, R, I ,L,V,K, 
2N,Y, YDOT, YL, YLDOT, III, COMPARO , COMPAR 1 , POT , V I NT , VI NC , MOMENTA , 
3DELTA0,DELTA1 ,XSEC0,XSEC1 , TOTAL, RA,ROS, AS, CHANGE, MANY, S,B,FNORM, BE 
4TTER , . 

DY( 1 ) = YC( 2 ) ' 

DY(2) = - (FK2(K) + V * ( EXPF ( -XC/A ) ) / ( XC/A ) ) *YC ( 1 ) 

RETURN 

END 



SUBROUTINE DER I VI ( XC , YC , DY ) 

DIMENSION E ( 100 ) ,FK( 100 ) ,FK2(100) , RN ( 100 ) , VAL JO ( 100 ) ,VALJ1( 100) 

1 ,VALJ2 ( 100) , VALNO(IOO), VALNl ( 100 ) ,VALN2 ( 100 ) , R(IOOO), Y(IOOO), 
2YDOT(1000), COMPAROI 100) , COMPAR 1 ( 1 00 ) , YL ( 1 000 ) , Y LDOT ( 1 000 ) , POT ( 1 
300) ,DELTAO( 100) ,DELTA1 (100) , XSECO ( 1 00 ) , XSEC 1 ( 100 ), TOTAL ( 100 ) ,RA(10 
40) ,FNORM( 100) , BETTER! 10,100) 

DIMENSION YC(30), DY(30) 

COMMON E,EINT,EINC, FK, FK2, RINT, RINC,RN, A, II, LL , NN , MODE, 

1 REND, VAL JO, VALJl ,VALJ2 , VALNO , VALN 1 , VALN2 , J J , M , NUMBER, R, I , L,V,K, 
2N,Y, YDOT, YL, YLDOT, III, COMPARO, COMPAR 1 , POT ,V I NT ,V I NC , MOMENTA , 
3DELTA0,DELTA1 ,XSEC0,XSEC1 , total, RA ,R0S, as, change, MANY, S,B,FN0RM, BE 
4TTER 

DY(1) = YC(2) 

DY(2) = - (FK2(K) + V * ( EXPF ( -XC/ A ) ) / ( XC/A ) -2./XC**2) * YC(1) 

RETURN 

END 
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SQUARE WELL POTENTIAL 



SUBROUTINE DERIV (XC»YC,DY) 

DIMENSION E ( 100 ) ,FK( 100) »FK2 ( 1 00 ) ♦ RN ( 1 00 ) »VALJ0( 100 ) »VALJ1 ( 100) 

1 »VALJ2( 100) ♦ VALNO(100)» VALNl ( 100 ) »VALN2 ( 100 ) ♦ R(IOOO), Y(IOOO)* 
2YDOT(1000), COMPARO( 100 ) , COMP AR 1 ( 1 00 ) ♦ YL (1000 ) »YLDOT ( 1000 ) »POT ( 1 
300) »DELTAO( 100) »DELTA1 ( 100 ) ,XSECO ( 100 ) »XSEC1 (100) ♦ TOTAL! 100) »RA(10 
40) »FNORM( 100 ) » BETTER (10.100) 

DIMENSION YC(30)» DY(30) 

COMMON E.EINT.EINC. FK, FK2, RINT. RING, RN » A, II. LL. NN . MODE. 
1REND.VALJ0.VALJ1.VALJ2.VALN0.VALN1.VALN2.JJ.M, NUMBER. R. I .L.V.K. 
2N.Y. YDOT. YL. YLDOT . III. COMPARO . COMPAR 1 . POT . V I NT . V I NC . MOMENTA . 
3 DELT AC. DELTA 1. XSECO. XSECl. TOTAL. RA.ROS, AS. CHANGE. MAN Y.S.B.FNORM. BE 
4TTER 

DY(1) =YC(2) 

IF(XC -A) 1.1.2 

1 DY(2) = -(FK2(K) + V) * YC ( 1 ) 

■ RETURN 

2 DY ( 2 ) =-( FK2 ( K) )* YC ( 1 ) 

RETURN 

END 



SUBROUTINE DER I VI ( XC . YC .DY ) 

DIMENSION E( 100 ) .FK( 100) .FK2 ( 100 ) . RN ( 100 ) .VALJ0( 100 ) .VALJK 100) 

1 .VALJ2( 100) . VALNO(IOO). VALNl ( 100 ) .VALN2 ( 100 ) . R(IOOO). Y(IOOO). 
2YDOT(1000). COMPARO! 100) . COMPAR 1 ( 1 00 ) . YL ( 1 000 ). Y LDOT ( 1000 ). POT ( 1 
300) .DELTAO! 100) .DELTAl ( 100 ) .XSECO ( 100) .XSECl (100) .TOTAL! 100) .RA( 10 
40) ,FNORM( 100 ) .BETTER ( 10.100 ) 

DIMENSION YC(30). DY(30) 

COMMON E.EINT.EINC. FK. FK2. RINT. RINC. RN. A. II. LL . NN. MODE. 

1 REND. VAL JO. VAL J 1.VALJ2.VALN0. VALNl. VALN2.JJ.M. NUMBER. R. I .L.V.K. 
2N.Y. YDOT. YL. YLDOT. III. COMPARO. COMPARl . POT . V I NT . V I NC . MOMEN T A , 
3DELTA0.DELTA1 .XSECO. XSECl. TOTAL. RA.ROS. AS.CHANGE. MAN Y.S.B.FNORM. BE 
4TTER 

DY(1) =YC(2) 

IF,(XC -A) 1.1.2 

1 DY(2) = -(FK2(K) + V - 2 . /XC**2 ) *YC ( 1 ) 

RETURN 

,2 DY(2) = -(FK2(K) - 2./ XC*#2 ) *YC(1) 

RETURN 

END 
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THE FOLLOWING DEFINITIONS FOR B MUST BE SUBSTITUTED IN READY2 FOR VARYING 



CATTERING LENGTH. THE DEFINITION PROVIDED IN READY2 IN THE LISTING IS FOR 



ARIATIONS IN EFFECTIVE RANGE. 



YUKAWA WELL 



B=ROS/ ( l.-1.369*SS+l .093*SS**2-1.127*SS**3+1 .005*SS**4-.324*SS**5 ) 



SQUARE WELL 



B = ROS/( l.-.5«-(S-l. ) + .486*(S-l. )**2-.476* (S-1. ) *^^3 + . 378* ( S-1 . )**4 



:-.158^f (S-1. )**5 ) 



exponential well 



B=ROS/( 1 .-.904*SS+.745*SS**2-.543*SS**3+.174*SS**4) 



GAUSSIAN WELL 

B=ROS/ ( l.-.641*SS+.568*SS**2-.721*SS**3+l. 172*SS**4-1.015*SS**5 ) 
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APPENDIX IV 



SAMPLE RESULTS 
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These results are for a yukawa potential with r 

I 

2.46 Fm and a varying in I 16% deviation from a * = 
s s 



* filled at 
,s 

23.677 FM. 
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¥ 



E F F E CT I V E ~“R A N G E / ‘SC ATT E R ITnrnrETrG T H 

-.8727457026E-01 
RADIUS POTENTIAL DEPTH 

.1086E-1 2 .5621E+.02 

ENERGY PHASE ANGLE CROSS SECTION, 



L = 0 



L=0 



.10 


.91 13E+00 


.3257E-22 


.20 


. 1037E+01 . 


.1933E-22 


.30 


. 1094E+01 


. 1 372E-22 


.40 


. 1 125E+01 


.1 062E-22 


.50 


. 1 145E+01 


.8649E-23 


.60 


. 1 158E+01 


.7290E-23 


.70 


. 1 166E+01 


.6295E-23 


.80 


. 1 172E+01 


.5535E-23 


.90 


. 1 176E+01 


.4936E-23 


1.00 


. 1 178E+01 


.4451 E-23 


1.10 


. 1 160E+01 


.4051 E-23 


1.20 


. 1 180E+01 


.3716E-23 


1.30 


. 1180E+01 


.3430E-23 


1 .40 


. 1 180E+01 


.3184E-23' 


1 .50 


. 1 179E+01 


.2970E-23 


1.60 


. 1 178E+01 


.2781 E-23 


1.70 


. 1176E+01 


.2S15E-23 


1 .80 


. 1 175E+01 


.2466E-23 


1 .90 


. 1 1 73E+01 


.2333E-23 


2.00 


. 1 172E+01 


.2213E-23 


2.10 


. 1 170E+01 


.2105E-23 


2.20 


. 1 168E + 01 - 


.2006E-23 


2.30 


. 1 166E+01 


.1915E-23 


2.40 


. 1 164E+01 


.18328-23 


2.50 


.1161E+01 


. 1755E-23 


2.60 


. 1 159E+01 


. 1685E-23 


ENERGY 


AJORMALIZED CROSS SECTION 


.10 


. 1 1366967E + C1 - 


.20 


■ . 1086931 1E+C1 


.30 


. 10658047E + C1 


.40 


.1054102 


OE + 01 


. 50 


. 10466564E + C1 


. 60 


. 1041 4947E + C1 


.70 


.10377CC 


2E + 01 


.80 


.1034789 


2L + 01 


.90 


.1032482 


1E + C1 


1.00 


. 10306062E + C1 


1.10 


. 1029C489E•^01 


1.20 


.1027733 


8E + 01 


1.30 


. 1026607 1E+C1 


.1 .40 


. 102563CCE + C1 


1.50 


.1024773 


6E + 01 


1 . 60 


. 1024016 1E+01 


1.70 


. 10233U05E+C1 


1 .80 


. 10227338E+C1 


1 .90 


.1022185 


3E+01 


2.00 


. 10216868E+01 


2.10 


. 1021 231 lE+Cl 


2.20 


. 102081 29E+01 


2.30 


. 10204272E+01 


2.40 


.1020070 


2E + 01 


2.50 


. 10197386E+C1 


2 ,J>.Q 


_ ^0J94296£ + 01 
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£FFECTIV£~RAKGE/SCATTERTI^G“ LENGTH 

-.91i43050217E-0K 
RADIUS POTENTIAL DEPTH 

.1G83E-12 .56M1E+02 

ENERGY PHASE ANGLE CROSS SECTION,' 



L=0 



L = 0 



.10 


.6913E+00 


.3155E-22 


.20 


. 1020E+01 


.1894E-22 


.30 


. 1079E+01 


.1351 E-22 


.40 


. 1 1 12E+01 


.1048E-22 


.50 


. 1 133E+01 


.8553E-23 


.60 


..1 146E+01 


.7218E-23 


.70 


. 1156E+01 


.6238E-23 


.80 


. 1 162E+01 


.5489E-23 


.90 


. 1 166E+01 


.4897E-23 


1.00 


. 1 169E+01 


.4419E-23 


1.10 


.1171E+01 


.40235-23 


1.20 


. 1172E+01 


.3691 E-23 


1.30 


. 1 172E+01 


.3403 E-23 


1 .40 


. 1 172E+01 


.3164E-23 


1 .50 


. 1172E+01 


.2952E-23 


1.60 


. 1171E+01 


.2765E-23 


1 .70 


. 1170E+01 


.26005-23 


1 .80 


. 1 168E+01 


.24535-23 


1.90 


.1167E+01 


.2321 £-23 


2.00 


. 1 165E+01 


.2202E-23 


2.10 


. 1 164E + 01 


.2094E-23 


2.20 


. 1 162E+01 


. 1995E-23 


2.30 


. 1 1 60E+01 


.19065-23 


2.40 


. 1 158E+01 


.1823E-23. 


2.50 


. 1 156E+01 


.1747 E-23 


2.60 


. 1 154E+01 


.1677 E-23 


ENERGY 


NORKALIZEC 


CROSS SECTION 


.10 


.11011434E+01 - 


.20 


. 10649739E+C1 


.30 


. 104936E6E 


+ 01 


.40 


. 10406591E+01 


. 50 


. 10350936E+01 


.60 


. 10312243E 


+ 01 


.70 


.10283743E+01 


.80 


. 1026 1847E+01 


.90 


. 10244473E+01 


1.00 


. 10230335E+01 


1.10 


. 1021 859CE+01 


1.20 


. 10208666E+01 


1.30 


. 102001 61E+01 


1.40 


. 10192782E+01 


1.50 


.10186313E+01 


1.60 


. 10180589E+01 


1.70 


. 10175484E+01 


1.80 


. 10170898E+C1 


1.90 


.10166753E 


+ 01 


2.00 


.10162984E+01 


•2.10 


. 10159539E+C1 


2.20 


. 10156377E+C1 


2.30 


.10153461E+01 


2.40 


. 10150762E+01 


2.50 


. 10148255E+01 


2-6Q_ 


.10J_4^919E + G1 
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EFFECTIVE RAKGE/SCJTTERTN^S LENGTH 

-.95586U3U09E-01 
RADIUS POTENTIAL DEPTH 

.1C80E-12 .5661E + 02 



ENERGY PHASE ANGLE CRCSS SECTICN, 




L = 0 ' 


L=0 


. TO 


.8719E+00 


.3056E-22 


.20 


. 1004E+01 


. 1855E-22 


.30 


. 1064E+01 


.1 329E-22 


.40 


. 1099E+01 


.1 034E-22' 


.50 


. 1121E+01 


.8457E-23 


.60 


.1135E+01 


.7 145E-23 


.70 


. 1 145E+01 


.6181 E-23 


.80 


. 1 152E+01 


.5443E-23 


.90 


. 1 157E+01 


.4859E-23 


1 .00 


. 1 161E+01 


.4386E-23 


1.10 


. 1163E+01 


.3994E-23 


1.20 


. 1 1o4E+01 


.3666E-23 


1.30 


. 1 165E+01 


.3386E-23 


1 .40 


. 1 165E+01 


.3144E-23 


1 .50 


. 1165E+01 


.2934E-23 


1.60 


. 1 164E+01 


.2749E-23 


1 .70 


. 1 163E+01 


.2585E-23, 


1.80 


. 1 162E+01 


.2439E-23I 


1 .90 


.1161 E+01 


.2308E-23: 


2.00 


. 1 1 59E+01 


.2190E-23 


2.10 


. 1 158E+01 


.20838-23 


2.20 


. 1 156E+01 


.1985E-23 


2.30 


. 1154E+01 


.1896E-23 


2.40 


. 1152E+01 


.1 814E-23 


2.50 


. 1 150 E+01 


.1733 E-23 


2.60 


. 1 148E+01 


. 1669E-23 


ENERGY 


KORiVArL IZ 


ED CROSS sect; 


. 10 


. 1066481 OE+01 


.20 


. 1043138CE+C1 


.30 


.103290 


3CE+C1 


.40 


. 1027 1472E+01 


.50 


. 10234527E+01 


.60 


. 10208768E + 01 


.70 


.10189756E+01 


.80 


.10175126E+01 


.90 


. 10163506E+01 


1.00 


.10154041E+C1 


1.10 


.10146173E+01 


1.20 


.10139521E+01 


1.30 


. 10133817E+01 


1.40 


. 10128867E + C1 


1 .50 


. 10124526E+01 


1.60 


. 10120684E+01 


1.70 


. 101 17257E+C1 


1.80 


.101141 


T8E+C1 


1 .90 


.1011 1394E+01 


2.00 


. 10108863E+01 


2.10 


. 10106549E + C1 


2.20 


. 10104425E+01 


2.30 


. 10102467E+01 


2.40 


. 10100654E+01 


2.50 


. 1009897CE + C1 


2.60 


.100974008+01^ 
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1 



EFFECTIVE R ANG E7 SC AT T ER I NG LENGTH 

-.9974236601 E-01 
RADIUS POTENTIAL DEPTF 

.1C77E-12 .5680E+C2 

ENERGY PHASE ANGLE CRCSS SECTION, 

L=0 L=0 



.10 


.8531E+00 


.2959E-22 


.20 


.9875E+00 


.1817E-22 


.30 


. 1050E+01 


.1 308E-22 


.40 


. 1086E+01 


.1021E-22 


.50 


. 1 1C9E+01 


.8360E-23 


.60 


. 1125E+01 


.7073E-23 


.70 


. 1135E+01 


.6124E-23 


.80 


. 1 143E+01 


.5396E-23 


.90 


. 1 148E+01 


.4820E-23 


1 .00 


. 1 152E+01 


.4352E-23 


1 . 10 


. 1 1 55E+01 


.3966E-23 


1.20 


.1156E+01 


.3641 E-23 


1.30 


. 1 157E+01 


.33638-23 


1 .40 


. 1 158E+01 


.3124E-23 


1 .50 


. 1 1 57E+01 


.29168-23 


1.60 


. 1 157E+01 


.2733E-23 


1.70 


. 1 156E+01 


.2570E-23 


1 .80 


. 1 155E+01 


.24258-23 


1 .90 


. 1 1 54E+01 


.2296E-23 


2.00 


.1153E+01 


.2 1788-23 


2.10 


. 1 152E+01 


.2072E-23 


2.20 


. 1150E+01 


.1975E-23 


2.30 


. 1148E+01 


.1886E-23 


2.40 


. 1 147E+01 


.1805 E-23 


2.50 


. 1 145E+01 


. 1730E-23 


2.60 


. 1 143E+01 


.1661E-23- 


ENERGY 


NORMALIZED 


CROSS SECTION 



. 10 


. 10327555E + 01 


.20 


. 1021 4671E + 01 


.30 


. 10164376E+C1 


.40 


.10135874E+01 


.50 


.101 17498E+01 


.60 


. 10104647E+C1 


.70 


. 10095142E+01 


.80 


. 100878 17E+01 


.90 


. 1008 1992E+01 


1.00 


. 10077243E + 01 


1.10 


. 10073292E+01 


1.20 


. 1006995CE+01 


1.30 


. 10067083E+01 


1 .40 


. 10064594E + C1 


1.50 


. 1006241 1E+01 


1.60 


.100604788+01 


1.70 


.100587538+01 


1.80 


. 100572C3E + C1 


1 .90 


.100558C2E+01 


2.00 


. 10054527E+01 


2.10 


. 10053363E + C1 


2.20 


. 10052293E + C1 


2.30 


. 10C513G7E+01 


2.40 


. 10050394E + 01 


2.50 


. 10049546E + 01 


. 2.60 


. 10048756E + C1 
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EFFECTIVE RANGE/SC ATT ERTTTG~rFl\GTH"" 

-.1C38982979E+00 
RADIUS POTENTIAL DEPTH 

.10748-12 .5700E+02 



ENERGY 


PHASE ANGLE CRCSS SECTICN, 




L=0 L 


= 0 


.10 


.8350E+00 


.2866 E-22 


.20 


.9717E+00 


. 1778E-22 


.30 


. 1036E+01 


. 1 287E-22 


.40 


. 1074E+01 


.1 007 E-22 


.50 


. 1098E+01 


.8263E-23 


.60 


. 1 1 14E+01 


.6999E-23 


.70 


. 1125E+01 


.6066E-23 


.80 


. 1 133E+01 


.5349E-23 


.90 


. 1 139E+01 


.4781 E-23 


1 .00 


. 1 143E+01 


.4319E-23 


1.10 


. 1 1 46E+01 


.3937E-23 


1.20 


. 1 148E+01 


.3615E-23 


1.30 


. 1 150E+01 


.3341 E-23 


1 .40 


. 1150 £+01 


.3104E-23 


1 .50 


. 1 1 SOE+01 


.2898E-23 


1.60 


. 1 150E+01 


.2716E-22 


1.70 


. 1 150E+01 


.2555E-23 


1.80 


. 1 149 £+01 


.2412E-23 


1 .90 


. 1 1 48E+01 


.2283E-23 


2.00 


. 1 147E+01 


.2 166 E-23 


2.10 


. 1 146E+01 


.20618-23 


2.20 


. 1 144E+01 


. 1965E-23 


2.30 


. 1143E+01 


.1 877E-23 


2.40 


. 1 141E+01 


.1796 E-23 


2.50 


. 1 140E+01 


.1721E-23 


2.60 


. 1138E+01 


.1653E-23-. 


ENERGY 


NORMALIZED 


CROSS SECTK 


.10 


. lOOOOOOOE+01 


.20 


. lOOOOOCCE+01 


.30 


. lOOOOOOCE+Cl 


.40 


.lOOOOOOOE 


+01 


.50 


. lOOOOOOOE+01 


.60 


. lOOOOOOCE+Ol 


.70 


. lOOOCOCCE+01 


.80 


. 1 OOOOOOOE+01 


.90 


. lOOOOOOOE+01 


1.00 


. lOOOOOCCE+Cl 


1.10 


. lOOOCCCGE+01 


.1.20 


. lOOOOOOOE+01 


1.30 


. lOOOOOOOE+Cl 


1 .40 


. lOOOOOCCE+01 


1 .50 


.lOOOOGCCE 


+ 01 


1.60 


. 1 OOOOOOOE+01 


1.70 


. lOOOOOOCE+01 


1.80 


. lOOOOOCCE+01 


1.90 


. lOOOOOCCE+01 


2.00 


. lOOOOOOOE+01 


2.10 


. lOOOOOOCE+Cl 


2.20 


. lOOOOOCCE+01 


2.30 


. 1 OOOCCCCE+01 


2.40 


. lOOOOOCOE+01 


2.50 


. lOOOOOOCE+01 


2.60 


. lOOOOOCCE+01 



EFFECTIVE' RANGE/"SC ATT ER I KG~rENGTH 

-.1C8C5U2298E+00 
RADIUS POTENTIAL DEPTH 

.1C71E-12 .5719E+C2 

ENERGY PHASE ANGLE CRCSS SECTICN, 



L=0 



L=0 



.10 


.8175E+00 


.2775E-22 


.20 


.9562E+00 


.1741E-22 


.30 


. 1022E+01 


.1266E-22 


.40 


. 1061E+01 


.9934E-23 


.50 


. 1086E+01 


.8166E-23 


.60 


. 1 103E+01 


.6926E-23 


.70 


.111 5E+01 


.6008E-23 


.80 


. 1 124E+01 


.5302E-23 


.90 


,11 30 E + 0 1 


.47418-23 


1 .00 


. 1135E+01 


.4286E-23 


1.10 


. 1 1 38E+01 


.3908E-23 


1.20 


. 1141E+01 


.3590E-22 


1.30 


. 1 142E+01 


.3318E-23 


1 .40 


. 1 143E+01 


.3084E-23' 


1 .50 


. 1 144E+01 


.2880E-23 


1 .60 


. 1 144E+01 


.2700E-23 


1 .70 


. 11U3E+01 


.25408-23 


1 .80 


. 1 143E+01 


.2398E-23 


1.90 


. 1 142E+01 


.2270E-23 


2.00 


.1141 £+01 


.21558-23 


2.10 


. 1 140E+01 


.2050E-23 


2.20 


. 1 1 39E+01 


. 1954E-23 


2.30 


. 1 1 37E+01 


.18678-23 


2.40 


. 1136E+01 


.1787E-22 


2.50 


. 1 134E+01 


.1713E-23 


2.60 


. 1 133E+01 


.1645E-23 • 


ENERGY 


NORMALIZED CROSS SECTK 


.10 


.96823662E + 00 


.20 


.97877150E+C0 


.30 


.98361 6C2E+CO 


.40 


.98640398E+00 


.50 


.9882181 1E+C0 


.60 


.98949447E+C0 


.70 


.99044257E+00 


.80 


.991 17554E+00 


.90 


.99175986E+00 


1.00 


.9922 


371 6E+C0 


1.10 


.99263482E+C0 


1.20 


.99297162E+00 


1.30 


.99326084E+00 


1 .40 


.9935 


1216E+C0 


1 .50 


.99373277E+00 


1.60 


.99392818E+00 


1.70 


.99410261E+00 


1.80 


.99425941E+00 


1 .90 


.9944C125E+G0 


2.00 


.99453026E+00 


2.10 


.99464820E+C0 


2.20 


.99475652E+C0 


2.30 


.9948564 1E+00 


2.40 


.99494888E+00 


2,50 


.99503479E+C0 


2.60 


.9951 1485E+00 
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EFFECTIVE RANGE /SC ATTER I NG LENGTH 

1122 1C16 18E + 00 
•RACIUS POTENTIAL DEPTH 

.1C68E-12 .5738.E±fi.2_ 

ENERGY PHASE ANGLE CRCSS SECTION, , 



L=0 



L = 0 



.10 


.80C5E+00 


.2686 E-22 


.20 


.94 12E+00 


. 1703E-22 


.30 


. 1009E+01 


. 1245E-22 


■ .40 


. 1049E+01 


.9797E-23 


. 50 


. 1075E+01 


.8068E-23 


. 60 


. 1093E+01 


.6852E-23 


.70 


. 1 105E+01 


.5950 £-23 


.80 


. 1 1 15E + 01 


.52548-23 


.'90 


. 1 122E+01 


.47C2E-23 


1 .00 


. 1 127E+01 


.4252E-23 


1.10 


. 1 130E+C1 


.3879E-23 


1.20 


. 1 133 E +01 


.3564E-23 


1.30 


. 1 135E+01 


.3296 E-23 


1 .40 


. 1 136E+01 


.3064E-23 


1 .50 


. 1 1 37E+01 


.2861 E-23 


1.60 


. 1 137E+01 


.2683E-23 


1.70 


. 1137E+01 


.2525E-23 


1 .80 


. 1 136E+01 


.2384E-23 


1 .90 


. 1 1 36E+01 


.2257E-23 


2.00 


. 1 135E+01 


.2 143E-23 


2.10 


. 1 154E+01 


.2039E-23 


2.20 


. 1 1 33E+01 


.1944E-23 


2.30 


. 1 1 32E+01 


.1 857E-23 


2.40 


. 1 130E+C1 


.1778E-23 


2.50 


. 1 129E+01 


.17048-23 


2.60 


. 1127E+01 


. 1636E-23 


ENERGY 


NORKALIZEC CROSS SECT 


.10 


.93747776E+00 


.20 


.95781212E+00 


.30 


.96730933E+C0 


.40 


.9728 1728E+00 


. 50 


.97641799E+00 


. 60 


.97895922E+C0 


.70 


.9808510 


3E+00 


.80 


.98231599E+00 


.90 


.9834853 


2E + C0 


1.00 


.9844414 


3E+C0 


1.10 


.98523866E+00 


1.20 


.9859143 


OE+00 


1.30 


■ .98649480E+Q0 


1.40 


.98699944E+C0 


1 .50 


.98744260E+00 


1.60 


.98783524E+00 


1.70 


.98818582E+00 


1.80 


.988501C4E+00 


1.90 


.98878622E+00 


2.00 


.9890456 


5E+00 


2.10 


.9892828 


5E + 00 


2.20 


.98950072E+GO 


2.30 


.98970165E+00 


2.40 


.98988768E+00 


2.50 


.9900605 


OE + CO 


2.60 


.99022157E+C0 
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EFFECTIVE RANGE/SCATTERING LENGTH 

-.11 63660957E+00 
RADIUS POTENTIAL DEPTH 



ENERGY 


PHASE ANGLE CRGSS SECTION, 




L=0 


L = 0 


.10 


.7842E+00 


.2601E-22 • 


.20 


.9264E+00 


.1667E-22 


.30 


.9957E+00 


.1224E-22 


.40 


. 1037E+01 


.9661E-23 


.50 


. 1064E+01 


.7971 E-23 


.60 


. 1082E+01 


.67788-23 


.70 


. 1096E+01 


.5892E-23 


.80 


. 1 1C6E+01 


.5207E-23 


.90 


.1113E+01 


.4662E-23 


1.00 


. 1 1 13E+01 


.4218E-23 


1 .10 


. 1 122E+01 


.3850E-23 


1.20 


. 1 125E+01 


.3539E-23 


1.50 


. 1 127E+01 


.3273E-23 


1 .40 


.1129E+01 ^ 


.3044E-23 


1.50 


. 1150 E +01 


.2843E-23 


1.60 


. 1 130E+01 


.26676-23 


1 .70 


. 1 130E+01 


.25106-23 


1.80 


. 1 130E+01 


.2370E-23 


1 .90 


. 1 1 30E+01 


.2244E-23 


2.00 


. 1 129E+01 


.2131E-23: 


2.10 


. 1 128E+01 


.2028 E-23 


2.20 


. 1127E+01 


. 1934E-23 


2.30 


. 1126E+01 


.1848 E-23 


2.40 


. 1 125E+01 


.1769E-23 ' 


2.50 


. 1 124E+01 


. 1696E-23 


2.60 


. 1 122E+01 


. 16286-23^ 


ENERGY 


NORHALIZEC CROSS SECTION 


■ .10 


.90772760E+C0 


.20 


.9371 4849E+00 


.30 


.951 1 01 69E+C0 


.40 


.95925670E+00 


.50 


.96461287E+C0 


;60 


.968404S9E+C0 


.70 


.971234 1 


8E+00' 


.80 


.9734287 


5E + 00 


.90 


.97518273E+00 


1.00 


.9766183 


2E + 00 


1.10 


.9778163 


4E+C0 


1.20 


.9788323 


2E+00 


1.30 


.97970572E+C0 


1 .40 


.98046532E+00 


1 .50 


.9811326 


3E+00 


1.60 


.9817240 


5E+00 


1.70 


.98225227E+C0 


1 .80 


.98272732E+00 


1 .90 


.98315716E+00 


2.00 


.93354827E+00 


2.10 


.98390590E+C0 


2.20 


.98423442E+00 


2.30 


.98453743E+C0 


2.40 


.98481798E+00 


2.50 


.98507864E+C0 


2.60 


...... .985321 5 aFj-QQ._ __ 
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EFFECTIVE RANGE/SCATTERING LENGTH 

-.12C5220256E+00 
RADIUS POTENTIAL DEPTH 

.1C62E-12 .5776E+C2 

ENERGY PHASE ANGLE CRCSS SECTION, 



L=0 



L = 0 



.10 


.7684E+00 


.2519E-22 


.20 


.9120E+00 


. 1630E-22 


.30 


.9828E+00 


.1203E-22 


.40 


. 1025E+01 


.9525E-23 


.50 


. 1053E+01 


.7873E-23 


.60 


. 1072E+01 


.67C4E-23 


.70 


. 1086E+01 


.5833E-23 


.80 


. 1097E+01 


.5159E-23 


.90 


. 1 1C4E+01 


.4622E-23 


1 .00 


. 1 nOE+01 


.4184E-23 i 


1 .10 


.111 4E+01 


.3320E-23 ■ 


1.20 


. 1 1 18E+01 


.3513E-23 


1 . 30 


. 1 120E+01 


.325GE-23 


1 .40 


. 1 122E+01 


.3023E-23 


1 .50 


. 1 123E+01 


.2825E-23 


1.60 


. 1124E+01 


.2650E-23 


1 .70 


. 1 124E+01 


.2495E-23 


1 .80 


. 1 124E+01 


.2356E-23 


1 .90 


. 1 124E+01 


.2231 E-23 


2.00 


. 1 123E+01 


.2 1 195-23 , 


2.10 


. 1 122E+01 


.2017E-23 


2.20 


. 1122E+01 


. 1923E-23 


2.30 


.1121 E+01 


.1 838E-23 : 


2.40 


. 1 1 19E+01 


.1760E-23 


2.50 


. 1 1 18 E+01 


.1637E-23 


2 .60 


. 1117E+01 


.1620E-23 - 


ENERGY 


NORKALIZEC CROSS SECTION 


. 10 


.8789833 


1E + 00 


.20 


.91680363E+C0 


.30 


.9350130 


CE+00 


.40 


.9457379 


5E+00 


.50 


• .9528 1525E + 00 


.60 


.957841 65E+C0 


.70 


.9616004 


5E+00 


.80 


.96452096E+00 


.90 


.9668581 9E+C0 


1.00 


.96377317E+C0 


1.10 


.97037256E+C0 


1.20 


.97 172985E+00 


1.30 


.97289732E+00 


1 .40 


.97391316E+C0 


1.50 


.9748C59 


2E+00 


1.60 


.975597415+00 


1.70 


.9763045 


3E + 00 


1 .80 


.9769406CE+C0 


1.90 


.9775 1627E+00 


2.00 


.97804014E+00 


2.10 


.97851924E+C0 


2.20 


.97895939E+00 


2.30 


.9793654 lE+00 


2.40 


.979741 3 


6E+00 


2.50 


.98009068E+CO 


2.60 , 


.98041627JL+C0 
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